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^Brad  H.  Nelson,  Ph.D. 

Co-PIs:  Rob  Holt,  Ph.D.,  John  Webb  Ph.D.,  Peter  Watson,  M.D. 

Title  of  Project:  Deciphering  the  Adaptive  Immune  Response  to  Ovarian  Cancer 

INTRODUCTION: 

Tumor-infiltrating  CD8+  T  cells  (CD8+  TIL)  are  strongly  associated  with  increased 
survival  in  ovarian  cancer.  However,  they  do  not  work  in  isolation.  Prior  to  our  DOD  Teal  Award, 
we  discovered  that  two  other  types  of  immune  cell  play  an  important  supportive  role:  B  cells  and 
helper  T  cells.  We  made  this  discovery  by  performing  a  systematic  analysis  of  immune  cells  in 
ovarian  cancer.  We  found  that  killer  T  cells  are  often  found  in  small  clusters  together  with  B  cells 
and  helper  T  cells.  Importantly,  we  found  that  patients  whose  tumors  have  these  combinations 
of  immune  cells  have  better  survival  rates  than  patients  whose  tumors  contain  killer  T  cells 
alone  (Figure  1).  This  tells  us  that  T  cells  and  B  cells  work  together  to  attack  tumors.  These 
findings  have  powerful  clinical  implications:  to  enhance  the  immune  response  to  ovarian  cancer, 
we  need  to  enhance  the  activity  of  all  three  types  of  immune  cell,  rather  than  killer  T  cells  alone. 

To  explain  these  observations,  we  hypothesized  that  B  cells  serve  as  “organizers”  that 
help  to  draw  T  cells  into  the  tumor.  In  addition,  we  reasoned  that  B  cells  might  present  tumor 
proteins  to  the  T  cells  to  facilitate  tumor  recognition.  We  further  hypothesized  that  helper  T  cells 
might  produce  cytokines  that  help  to  stimulate  the  killer  T  cells.  To  test  these  hypotheses,  in  our 
Teal  Award  we  proposed  to  determine  which  tumor  proteins  (antigens)  are  recognized  by  B 
cells  and  helper  T  cells  in  high-grade  serous  ovarian  cancer  (HGSC).  By  identifying  these 
antigens,  we  would  create  new  molecular  tools  to  elucidate  how  the  immune  system  recognizes 
and  attacks  ovarian  cancer.  The  study  had  four  tasks: 

Task  1 .  To  identify  tumor  antigens  recognized  by  CD20+  TIL. 

Task  2.  To  identify  tumor  antigens  recognized  by  CD4+FoxP3+  TIL. 

Task  3.  To  determine  whether  tumor-infiltrating  B  cells  and  T  cells  recognize  the  same 
antigens. 

Task  4.  To  assess  the  functional  phenotype  of  antigen-specific  CD4+FoxP3+  TIL. 

Significance:  The  immune  system  has  a  profound  influence  on  survival  from  ovarian 
cancer.  With  better  understanding  of  the  immune  response,  it  will  be  possible  to  design  new 
treatments  such  as  vaccines  and  adoptive  T  cell  therapies  that  enhance  tumor  immunity  and 
increase  patient  survival.  We  envisioned  our  work  would  lead  to  a  major  re-think  about  cancer 
vaccines:  instead  of  simply  trying  to  activate  killer  T  cells,  we  also  need  to  find  effective  ways  to 
activate  their  team  mates,  the  B  cells  and  helper  T  cells. 

BODY: 

Task  1.  To  identify  tumor  antigens  recognized  by  CD20+  TIL. 

In  this  task,  we  proposed  to  identify  the  antigens  recognized  by  the  3  most  abundant 
CD20+  TIL  clones  from  each  of  3  ovarian  cancer  patients.  To  accomplish  this,  we  proposed  to 
clone  immunoglobulin  G  (IgG)  molecules  from  individual  CD20+  TIL.  These  would  be  used  to 
identify  the  corresponding  antigens  using  three  different  approaches:  candidate  antigen  assays, 
cDNA  library  screening,  and  mass  spectrometry. 
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Progress: 

1.1  Generating  IpG  and  TCR  profiles  from  ovarian  tumor  specimens 

To  identify  the  3  most  abundant  IgG  (for  Task  1 )  and  TCR  clones  (for  Task  2)  in  each 
patient’s  tumor,  we  used  deep  sequencing  methods  developed  in  the  Holt  lab.  Our  approach  for 
deep  sequence  analysis  of  T  cell  receptor  (TCR)  beta  chain  diversity  has  been  previously 
described  (1, 2),  and  we  made  adaptations  to  enable  B  cell  receptor  (BCR)  heavy  chain  profiling 
and  sample  multiplexing,  as  described  below.  We  isolated  total  RNAfrom  bulk  ascites  from  the 
three  HGSC  ovarian  cancer  study  subjects  at  three  different  timepoints  (the  time  of  primary 
disease,  at  first  recurrence,  and  at  second  recurrence).  First  strand  cDNA  was  synthesized 
from  total  RNA  using  either  TCR  beta  chain  or  BCR  heavy  chain  gene-specific  primers,  and  this 
cDNA  was  then  used  as  template  to  amplify  the  CDR3  sequences  from  each  receptor 
population  by  PCR.  The  CDR3  region  is  the  site  of  VDJ  recombination  and  is  the  most  highly 
variable  and  thus  most  informative  region  for  the  purpose  of  profiling  repertoire  diversity.  The 
PCR  primers  were  tagged  with  unique,  6  base  pair,  error  tolerant  barcode  sequences.  A 
different  barcode  was  used  for  each  individual  sample.  The  CDR3  amplicons  were  then 
subjected  to  several  additional  rounds  of  nested  PCR  using  primers  tailed  with  lllumina  adapter 
sequences  necessary  for  annealing  to  the  lllumina  flow  cell  surface  for  massively  parallel 
sequencing.  The  samples  were  pooled  and  sequenced  using  a  single  flow  cell  and  the  lllumina 
miSeq  platform.  Sequence  data  was  then  de-convoluted  based  on  the  barcode  identifiers,  and 
the  TCR  and  BCR  CDR3  sequences  were  filtered  to  remove  low  quality  reads,  then  clustered 
using  standard  bioinformatic  methods.  Clustering  was  done  to  determine  the  number  of  distinct 
CDR3  sequences  present  in  each  sample,  which  reflects  the  diversity  of  distinct  T-cell  and  B- 
cell  clonotypes  with  the  original  ascites  specimens.  Prior  to  reverse  transcription  and  PCR  each 
sample  was  spiked  with  RNA  from  a  single,  known  clonotype  in  order  to  monitor  PCR  and 
sequence  error  rates  and  the  efficiency  of  sequence  recovery. 

We  filtered  and  analyzed  the  sequence  data  derived  from  specimens  taken  from  each  of 
the  three  timepoints  for  one  of  the  three  study  subjects.  We  obtained  over  30,000  total  TCR 
beta  chain  CDR3  sequences  from  each  of  these  three  samples,  and  over  15,000  total  BCR 
heavy  chain  CDR3  sequences  from  each  of  these  same  three  samples.  These  collapsed,  upon 
clustering,  into  a  much  smaller  number  of  distinct  sequences.  Specifically  there  were  1 ,598 
distinct  TCR  sequences  and  688  distinct  BCR  sequences  observed  per  sample.  However,  as 
expected,  the  majority  of  these  unique  sequences  were  “singletons”,  representing  rare 
clonotypes,  with  the  bulk  of  the  accumulated  sequence  data  from  each  sample  represented 
approximately  30  clonotypes  per  sample.  Thus  the  ascites  repertories  are  polyclonal,  but  not  to 
the  same  extent  as  peripheral  blood,  where  previously  we  have  been  able  to  detect  over  a 
million  distinct  sequences  from  a  single  sample  by  deep  sequencing.  Interestingly,  18  of  the 
BCR  clonotypes  were  detectable  at  all  three  timepoints,  albeit  with  varying  abundance,  but  the 
TCR  clonotypes  were  much  more  sample  specific,  with  only  one  TCR  clonotype  being  present 
in  all  three  samples. 

1.2  Amplification  of  matched  IpG  heavy  and  light  chains  by  single-cell  RT-PCR 

Single  cell  sorting  and  PCR  amplification  of  matched  immunoglobulin  (Ig)  heavy  and 
light  chains  from  tumor-infiltrating  B  lymphocytes  was  accomplished  using  methods  similar  to 
Tiller  et  al.  (3).  Single  B  cells  were  isolated  based  on  their  surface  expression  of  CD19,  CD20 
and  IgG  using  a  BD  Influx  cell  sorter.  Individual  cells  were  sorted  directly  into  single  wells  of  96- 
well  PCR  plates.  Reverse  transcription  was  accomplished  directly  within  the  sorting  plate,  which 
yielded  cDNA  corresponding  to  the  original  single  cell.  This  cDNA  was  then  split  equally 
between  four  separate  PCR  reactions,  one  each  for  the  variable  portions  of  Ig  heavy,  kappa, 
and  lambda,  and  one  for  the  control  housekeeping  gene  GAPDH.  The  variable  portions  of  the 
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immunoglobulin  genes  were  then  amplified  by  sequential  nested  multiplex  PCR  using  primer 
sets  known  to  amplify  all  human  variable  and  joining  gene  segments  and  containing  restriction 
enzyme  cut  sites  to  facilitate  downstream  molecular  cloning.  We  successfully  amplified  matched 
pairs  of  Ig  genes  from  B  cells  isolated  from  the  tumors  of  three  high-grade  serous  ovarian 
cancer  patients.  Roughly  70%  of  all  sorting  wells  yield  amplification  of  at  least  one  B  cell- 
specific  gene  product,  while  about  10-25%  yield  matched  pairs  of  heavy  and  light  chains.  In 
total,  we  successfully  amplified  over  30  pairs  of  Ig  genes  from  the  3  patients. 

1.3  Cloning  of  CDR3  regions  into  IpG  expression  vectors 

Molecular  cloning  of  the  PCR  amplified  variable  portions  of  matched  heavy  and  light 
chain  genes  derived  from  single  sorted  tumor  infiltrating  B  cells  was  accomplished  by  standard 
techniques.  All  amplified  pairs  were  subsequently  cloned  into  appropriate  expression  vectors 
containing  the  signal  sequence  and  constant  portions  of  the  corresponding  germline 
immunoglobulin  chain  (heavy,  kappa,  lambda).  Once  inserted  into  these  expression  vectors,  the 
DNA  sequence  of  each  immunoglobulin  chain  was  determined  by  Sanger  sequencing. 
Sequences  for  each  of  the  pairs  of  heavy  and  light  chains  were  obtained.  We  evaluated  the 
prevalence  of  each  immunoglobulin  heavy  sequence  by  sequencing  samples  prepared  by 
conventional  RT-PCR  of  bulk  tumor  preparations. 

1.4  Preparation  of  recombinant  IpG  (rIpG) 

To  produce  recombinant  IgG,  two  expression  vectors  encoding  cloned  matched  heavy 
and  light  chain  variable  portions  were  co-transfected  into  log-phase  293T  cells.  This  transient 
transfection  was  accomplished  by  a  calcium  phosphate  precipitation  method.  After  seven  days 
in  culture,  supernatants  were  harvested  and  assessed  for  the  presence  of  rIgG  by  human  IgG- 
specific  western  blot.  Fully  assembled  IgG  molecules  (-150  kDa)  were  observed  in  the 
supernatant  of  80-90%  of  culture  transfected  with  different  pairs  of  expression  vectors  under 
non-reducing  conditions.  To  isolate  rIgG  from  the  culture  supernatant,  it  was  concentrated  and 
run  over  protein  G  sepharose  columns.  Captured  IgG  molecules  were  eluted  at  low  pH  and  the 
resulting  eluate  was  neutralized  and  desalted.  Further  concentration  and  assessment  of  purity 
was  accomplished  by  BCA  assay  and  additional  western  blot  analysis.  Over  30  recombinant 
antibodies  were  successfully  expressed  in  vitro,  for  an  overall  success  rate  of  80%. 

1.5  Discovery  of  plasma  cells  as  the  major  source  of  IpG  in  HGSC 

Unexpectedly,  none  of  the  IgG  sequences  from  the  above  CD20+  B  cells  aligned  to  the 
dominant  clonotypes  determined  by  sequencing  of  bulk  tumor  tissue.  Furthermore,  on  screening 
recombinant  antibodies  against  HGSC  cell  lines,  we  obtained  no  evidence  of  tumor  reactivity. 
This  puzzling  result  led  us  to  hypothesize  that  the  B  cell  component  of  TIL  might  be  more 
heterogeneous  than  previously  appreciated.  We  therefore  turned  to  immunohistochemistry  and 
flow  cytometry  to  further  investigate  B-lineage  cells  in  HGSC. 

To  this  end,  we  performed  multicolor  immunohistochemistry  on  whole  tumor  sections 
from  21  HGSC  patients.  We  found  that  CD20+  B  cells  often  form  dense  aggregates  with  other 
lymphocytes  within  and  around  HGSC  tumors  (Figures  2  and  3).  These  aggregates  resemble 
lymph  nodes,  yet  are  embedded  within  tumor  tissue.  Others  have  named  these  “ectopic 
lymphoid  follicles”  or  “tertiary  lymphoid  structures”  (TLS);  we  will  use  the  latter  designation.  TLS 
appear  to  have  all  the  components  of  normal  lymph  nodes.  They  have  T  cell  zones  containing 
CD4+  and  CD8+T  cells;  activated  dendritic  cells  expressing  CD208;  and  high  endothelial 
venules  (HEVs).  They  also  contain  B  cell-rich  areas  (called  follicles),  which  contain 
interdigitating  networks  of  follicular  DCs,  follicular  helper  T  cells  and  germinal  center  B  cells 
(Figure  3).  We  found  the  latter  cells  express  the  transcription  factor  Bcl-6  and  activation- 
induced  cytidine  deaminase  (AID),  an  enzyme  involved  in  class-switch  recombination  and 
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somatic  hypermutation.  Thus,  TLS  resemble  lymph  nodes  in  their  organization,  cellular 
composition,  and  ability  to  support  immune  reactions. 

By  flow  cytometry,  we  found  that  the  vast  majority  of  tumor-infiltrating  CD20+  B  cells  are 
lgG+  memory  B  cells  (Figure  4).  Most  of  these  memory  B  cells  appear  to  be  activated  (i.e.,  they 
express  FAS)  and  express  IgG.  However,  we  also  discovered  that  9/10  cases  contained  an 
appreciable  number  of  early  plasma  cells  (CD38'^'®^,  CD138'°",  IgD-)  (Figure  4).  The  proportion 
of  these  cells  ranged  from  <5%  to  >80%  of  GDI  9+  cells.  Importantly,  these  plasma  cells  do  not 
express  CD20  and  therefore  have  been  missed  in  previous  immunohistochemical  and  sorting 
experiments  (Figure  5).  Plasma  cells  are  known  to  express  large  amounts  of  IgG-encoding 
mRNA  and  protein.  In  fact,  they  are  often  referred  to  ask  “antibody  factories”.  Given  this,  we 
hypothesized  that  IgG  sequencing  from  bulk  tumor  preparations  might  largely  reflect  the  IgG 
repertoire  of  plasma  cells.  To  address  this,  we  used  conventional  FACS  to  isolate  CD20+  TIL 
and  CD38'^'®'^  plasma  cells  from  tumor  samples.  We  then  sequenced  the  IgG  repertoire  in  each 
sample  and  compared  these  sequences  to  those  found  in  bulk  tumor  preparations  from  the 
same  tumor  (Figure  6).  In  3/3  patients,  we  found  that  both  CD20+  TIL  and  CD38'^'®^  plasma 
cells  were  clonally  expanded  to  some  extent.  We  also  found  examples  of  shared  IgG  sequences 
between  these  two  populations,  suggesting  that  some  CD20+  TIL  are  products  of  the  same 
immune  reactions  that  give  rise  to  plasma  cells.  Importantly,  we  found  that  IgG  sequences 
derived  from  plasma  cells  were  highly  over-represented  in  sequencing  data  from  bulk  tumor 
(Figure  6).  Thus,  we  concluded  that  plasma  cells  in  HGSC  are  the  major  source  of  clonally 
expanded  IgG  sequence  and  therefore  are  the  more  likely  B  lineage  cell  type  to  mediate  anti¬ 
tumor  immune  reactions. 

A  manuscript  describing  the  above  findings  is  under  review  at  the  AACR  journal 
Clinical  Cancer  Research  (Appendix  A). 

1.6  Cloning  and  production  of  recombinant  IpG  from  tumor-derived  plasma  cells 

Having  identified  plasma  cells  as  the  major  source  clonally  expanded  B  cell  populations 
in  HGSC,  we  investigated  the  possibility  that  these  cells  are  making  antibodies  to  tumor- 
associated  antigens.  As  per  the  original  proposal,  we  used  single-cell  sorting  followed  by  RT- 
PCR  to  molecularly  clone  matched  IgG  heavy  and  light  chains  from  these  cells.  We  isolated  and 
PCR-cloned  15  matched  IgG  heavy  and  light  chain  pairs  from  plasma  cells  from  an  HGSC 
patient  (Figure  6).  As  expected,  some  of  these  antibodies  correspond  to  clonally  expanded 
populations  dominating  the  IgG  repertoire  in  both  the  plasma  cell  component  and  the  bulk  tumor 
(Figures  6  and  7).  We  expressed  these  cloned  heavy  and  light  chains  as  recombinant 
antibodies  (n=15)  (Figure  6). 

1.7  Current  work:  Identifying  the  antipens  recognized  by  plasma  cell-derived  IpG 

With  15  recombinant  antibodies  in  hand,  and  an  efficient  pipeline  for  producing  more,  we 
are  in  a  strong  postion  to  identify  the  underlying  tumor  antigens.  Unfortunately,  we  did  not  reach 
this  milestone  during  the  funding  period  of  this  grant.  However,  we  are  actively  pursuing  this 
goal  using  two  new  sources  of  leveraged  funding: 

(a)  With  a  new  operating  grant  from  the  Canadian  Institutes  of  Health  Research  (CIHR) 
(see  Leveraged  Funding),  we  are  using  an  innovative  protein  array-based  approach  to  antigen 
identification.  Specifically,  our  approach  is  to  screen  microarrays  spotted  with  approximately 
9,000  cancer-  and  non-cancer-associated  proteins  (ProtoArray,  Life  Technologies),  as  others 
have  done  in  other  cancers  and  autoimmune  conditions  (90-92).  Importantly,  the  ProtoArray 
includes  cell  surface  proteins  (e.g.,  ErbB2)  that  have  been  shown  to  be  recognized  by  serum 
antibodies  from  ovarian  cancer  patients  (93).  Protein  arrays  will  be  pre-blocked  and  probed  with 
diluted  rIgG.  Bound  rIgG  will  be  detected  with  AlexaFluor647-conjugated  anti-human  IgG 
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secondary  antibody.  Arrays  will  be  scanned  using  our  LI-COR  Odyssey,  and  immunoreactive 
antigens  will  be  identified  using  ProtoArray  Prospector  v5.2.  Key  results  will  be  validated  by 
ELISA  using  recombinant  antigens  (Bacpower™  service;  Genscript). 

(b)  With  a  new  pilot  grant  from  Canada’s  Networks  of  Centres  of  Excellence  (NCE) 
program,  we  are  collaborating  with  Dr.  John  Babcook  at  the  Centre  for  Drug  Research  and 
Development  (CDRD)  at  the  University  of  British  Columbia  in  Vancouver.  Dr.  Babcook  is  an 
international  expert  in  antibody  production  and  engineering.  He  is  using  mass  spectrometry- 
based  approaches  (proteomics)  to  identify  the  antigens  recognized  by  the  15  plasma  cell- 
derived  antibodies  described  above. 

Thus,  we  remain  committed  to  the  goal  of  identifying  the  tumor  antigens  recognized  by 
tumor-associated  plasma  cells  in  HGSC.  Based  on  the  strength  of  the  data  generated  by  our 
DOD  Teal  Award,  we  have  been  successful  at  obtaining  continued  funding  for  this  challenging 
but  important  objective. 

Task  2.  To  identify  tumor  antigens  recognized  by  CD4+FoxP3+  TIL. 

We  proposed  to  clone  TCR  molecules  from  CD4+FoxP3+  tumor-infiltrating  T  cells  from 
the  3  ovarian  cancer  patients  described  above.  Cloned  TCR’s  were  to  be  expressed  in  a 
hybridoma  cell  line,  which  in  turn  would  be  used  to  identify  the  underlying  antigens  using  a 
candidate  approach  and/or  cDNA  library  screening. 

Progress: 

2.1  Discovery  of  extensive  clonal  diversity  among  tumor-infiltrating  T  cells 

For  Task  2,  we  originally  proposed  to  choose  from  ascites  from  each  of  3  patients  with 
EOC  (i.e.,  high-grade  serous  ovarian  cancer),  the  3  most  abundant  and  persistent  CD4+FoxP3+ 
TIL  clones  as  identified  by  TCR-seq  analysis,  and  proceed  to  identify  their  cognate  antigens. 
TCR-seq  analysis  was  completed  but  gave  the  surprising  result  that  ascites  T-cells  were  highly 
heterogeneous  and  we  did  not  find  any  examples  of  abundant,  persistent  T-cell  clones  that  were 
present  in  all  three  samples  from  a  given  subject.  If  we  had  observed  a  small  number  of 
dominant  clontypes  it  would  have  made  sense  to  proceed  with  conventional  cloning  and  antigen 
identification  as  described  in  our  original  proposal.  Now,  however,  we  are  faced  with  high  clonal 
diversity  of  tumor  associated  T  cells  and  the  question  of  what  antigens  these  populations  of  T 
cells  are  recognizing  remains.  This  prompted  us  to  refocus  our  approach,  and  work  towards 
developing  a  more  powerful,  high-throughput  and  unbiased  method  for  T  cell  antigen  discovery. 

2.2  Development  of  a  novel  high  throughput  screen  to  identify  antipens  recognized  by  T  cells 

As  we  have  discussed  in  recent  review  articles,  it  is  well  recognized  that  TCR-seq  has 
opened  a  new  window  on  T-cell  repertoires,  but  T-cell  antigen  profiling  capabilities  have  not 
kept  pace.  The  difficulties  associated  with  scaling  T-cell  antigen  profiling  stem  from  the 
enormous  complexity  of  the  pMHC/TCR  interaction.  This  complexity  is  generated  by  differential 
peptide  processing  in  antigen-presenting  cells,  MHC  allelic  variation,  somatic  recombination  of 
the  TCR  loci  in  T-cells,  and  the  requirement  for  co-receptors  (CDS  or  CD4)  in  these  interactions. 

Using  the  model  mouse  T-cell  antigen,  chicken  ovalbumin,  and  a  cytolytic  T-cell  (CTL) 
receptor  specific  for  this  antigen  (derived  from  the  TCR-transgenic  mouse  model,  OT-I)  we  have 
developed  and  validated  a  screening  protocol  and  data  analysis  pipeline  for  large  scale 
detection  of  CTL  epitopes.  The  primary  screening  consists  of  co-incubating  T-cells  of  interest 
and  target  cells.  These  target  cells  have  been  modified  to  carry  short  peptide  coding  sequences 
linked  to  a  FRET  reporter  construct  composed  of  CFP  and  YFP  monomers  separated  by  a 
granzyme-B-cleavable  linker.  When  cognate  peptides  are  recognized  by  T-cell  receptors  of 
interest,  granzyme-B  is  transferred  into  the  target  cells  carrying  the  immunogenic  peptide.  This 
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results  in  cleavage  of  the  FRET  reporter,  which  produces  a  shift  in  the  optical  properties  of  the 
cell,  and  allows  us  to  select  for  these  cells  using  fluorescence-activated  cell  sorting  (FACS). 
Collected  cells  are  then  lysed  and  the  lentiviral  transgene  is  amplified  for  sequencing  by  next- 
generation  methods  to  reveal  the  identity  of  peptide-coding  sequences  responsible  for  eliciting 
CTL-reactivity. 

The  specificity  of  the  approach  was  validated  by  creating  a  positive-control  cell  line 
expressing  the  correct  ovalbumin  epitope  (IDS.minOva)  and  a  decoy  cell  line  coding  for  a 
scrambled  version  of  the  epitope  (IDS.minOScr).  When  mixed  at  a  1:1  ratio  and  subjected  to 
OT-I  CTL  interrogation,  nearly  complete  enrichment  of  the  correct  sequence  was  observed  in 
the  genomic  DNA  of  collected  cells  (Figure  8A).  To  assess  the  sensitivity  of  the  method,  a  high- 
diversity  library  of  completely  random  peptide-coding  sequence  was  constructed.  The  positive 
control  line  (IDS.minOva)  was  spiked  in  to  the  random  minigene  encoding  cell  line  at  varying 
abundances  and  was  similarly  subjected  to  screening  with  OT-I  CTL.  With  minimal  optimization 
and  only  a  single-round  of  screening  (no  additional  panning  rounds),  Ova-coding  sequences 
were  detected  above  background  at  abundances  as  low  as  100  ppm  (Figure  8B). 

2.3  Current  work:  identifying  the  antipens  recognized  by  TIL  in  HGSC 

Based  on  the  strong  data  generated  by  our  DOD  Teal  Award,  we  were  successful  in 
obtaining  new  funding  from  Genome  Canada,  Genome  British  Columbia,  and  BC  Cancer 
Foundation  for  our  T  cell  antigen  identification  efforts  (see  Leveraged  Funding).  Our  next  steps 
are  to  apply  the  above  approach  to  detecting  cognate  antigens  from  tumor  cDNA  from  the 
ovarian  cancer  subjects  studied  here  as  well  as  to  characterize  TCR  cross-reactivity  of  tumor- 
reactive  T-cells.  We  are  starting  by  evaluating  CDS'"  cytolytic  T  cells  from  patient  ascites  or  bulk 
tumor.  We  will  isolate  CD45^from  samples  and  expand  these  isolated  lymphocytes  ex  vivo. 

From  CD45'  tumor  cells,  cDNA  fragment  libraries  will  be  generated  and  used  to  create  a  library 
of  cDNA-encoding  lentivirus,  which  will  be  transduced  into  autologous  B-cells  (from  patient 
blood  samples).  Transduced  B-cells  will  then  be  subject  to  the  screening  and  sequencing 
approach  that  has  been  validated  above.  Peptides  encoded  by  positive  hits  will  then  be  used  to 
selectively  enrich  tumor-infiltrating  or  tumor-associated  lymphocyte  populations  for  use  as 
increased-potency  adoptive  cell  therapeutics.  T-cell  cross-reactivity  is  another  important 
consideration  that  we  will  examine  by  using  the  above  method  to  screen  patient-derived  T-cell 
clonal  lines  of  interest  against  completely  random  peptide  coding  sequence.  Though  cross¬ 
reactivity  is  typically  not  observed  in  healthy  individuals  due  to  natural  mechanisms  of  central 
and  peripheral  tolerance,  mounting  evidence  has  indicated  that  each  T-cell  receptor  is  capable 
of  recognizing  a  vast  array  of  unique  sequences  in  complex  with  their  complementary  MHC 
molecule(s).  Using  random  sequence  that  is  not  restricted  to  naturally  occurring  coding 
sequence,  we  will  assess  the  level  of  cross-reactivity  observed  in  highly  tumor-reactive  T-cells 
and  work  to  develop  a  predictive  metric  that  will  provide  insight  into  the  probability  of  adverse 
events  resulting  from  administration  of  highly  cross-reactive  T-cells  versus  relatively  lower 
cross-reactive  T-cells. 

Task  3.  To  determine  whether  tumor-infiltratincf  B  cells  and  T  cells  recognize  the  same 
antigens. 

The  goal  of  this  task  was  to  test  the  hypothesis  that  CD20+  TIL  can  serve  as  antigen 
presenting  cells  (ARC)  in  the  tumor  microenvironment. 

Progress: 

As  discussed  in  Task  1  (above),  with  this  DOD  Teal  Award  we  learned  that  CD20+  TIL 
are  highly  polyclonal  (compared  to  their  plasma  cell  counterparts).  This  makes  it  far  less  likely 
that  they  serve  as  ARC  in  the  tumor  microenvironment.  Instead,  we  are  exploring  the  potential 
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role  of  plasma  cells  in  this  process  (as  described  in  Task  1).  As  for  CD20+  TIL,  we  are  now 
investigating  the  alternative  hypothesis  that  they  serve  as  “organizers”  of  lymphocytic  infiltrates 
by  releasing  immune  stimulatory  cytokines  and  chemokines  in  the  tumor  bed.  With  new, 
leveraged  funding  from  the  Canadian  Institutes  of  Health  Research  (CIHR)  (see  Leveraged 
Funding),  we  are  using  an  innovative  platform  (Nanostring)  to  perform  gene  expression 
profiling  of  HGSC  tumors  that  contain  either  (a)  CD8+  TIL  alone,  or  (b)  CD8+  and  CD20+  TIL 
together.  From  this  work,  we  intend  to  identify  immune  stimulatory  factors  released  by  CD20+ 
TIL  that  enhance  the  anti-tumor  effects  of  CD8+  TIL. 

Task  4.  To  assess  the  functional  phenotype  of  antigen-specific  CD4+FoxP3+  TIL. 

CD4+FoxP3+  TIL  show  great  functional  heterogeneity  in  EOC,  which  makes  it  difficult  to 
draw  definitive  conclusions  about  their  role  in  the  tumor  environment.  With  knowledge  of  their 
cognate  antigens,  we  will  be  able  to  clarify  this  issue  by  assessing  the  functional  phenotype  of 
individual  CD4+FoxP3+  T  cell  clones  as  opposed  to  bulk  cell  preparations.  We  proposed  to  do 
this  by  constructing  MHC  class  II  tetramers,  which  bind  specifically  to  CD4+  T  cells  expressing 
TCRs  relevant  to  a  particular  antigen. 

Progress: 

In  flow  cytometry  experiments,  we  discovered  that  CD4+FoxP3+  TIL  comprised  two 
distinct  populations  that  could  be  distinguished  by  their  expression  of  CD25.  In  fact,  the 
entire  CD4+  TIL  compartment  could  be  subdivided  into  four  subpopulations  based  on 
expression  of  FoxP3  and  CD25  (Figure  9).  We  performed  additional  flow  cytometry 
experiments  to  assess  the  activation  status  and  cytokine  production  profiles  of  these  four 
CD4+  TIL  subsets  (n  =  6  cases).  CD25+FoxP3+  TIL  had  a  classic  regulatory  T  cell 
phenotype  (Figure  9).  Unexpectedly,  we  discovered  that  the  CD25+FoxP3-  TIL 
subpopulation  had  a  previously  undescribed  phenotype.  On  the  one  hand,  they  had  an 
activated  phenotype  as  evidenced  by  high  CD69  and  low  CCR7  expression.  In  comparison 
to  Tregs,  CD25+FoxP3-  TIL  expressed  similar  levels  of  GITR;  however,  they  expressed 
lower  levels  of  CTLA-4  and  0X40  and  were  negative  for  Helios  (Figure  9).  These  data 
suggested  that  CD25+FoxP3-  cells  might  be  Thi  cells,  which  are  widely  reported  among 
TIL.  However,  CD25-+FoxP3-  T  cells  failed  to  produce  any  of  the  hallmark  Thi  cytokines 
IFN-y,  TNF-a  or  IL-2  after  in  vitro  stimulation  with  PMA  and  ionomycin  (Figure  9C).  Indeed, 

Thi  cytokines  were  only  produced  by  the  CD25-FoxP3-  subset.  None  of  the  CD4+  TIL 
subsets  produced  IL-4  or  IL-17A  (not  shown). 

Given  that  CD25+FoxP3-  T  cells  did  not  express  canonical  Th  cytokines,  we 
investigated  the  possibility  that  they  might  represent  other  less  common  CD4+  T  cell 
phenotypes.  CD4+  cytolytic  T  cells  have  been  described  in  cancer;  however,  similar  to  other 
non-Treg  cells,  very  few  CD25+FoxP3-  T  cells  expressed  the  cytolytic  markers  TIA-1, 
granzyme  B  or  perforin  (not  shown).  Furthermore,  CD25+FoxP3-  TIL  did  not  express 
CXCR5,  a  marker  of  T  follicular  helper  cells  (not  shown). 

Based  on  their  lack  of  discernible  functional  attributes,  we  hypothesized  that 
CD25+FoxP3-TIL  might  be  in  a  suppressed  or  exhausted  state.  In  accord  with  this,  we 
found  that  CD25+FoxP3-  T  cells  expressed  very  high  levels  of  the  exhaustion  marker  PD-1 
(Figure  9D).  Indeed,  the  level  of  PD-1  expressed  by  CD25+FoxP3-  TIL  was  on  average  3.1- 
fold  higher  than  Tregs,  and  6.6-fold  higher  than  CD25-FoxP3-  cells.  In  addition, 

CD25+FoxP3-  TIL  expressed  the  exhaustion  markers  LAG-3  (not  shown)  and  TIM-3,  the 
latter  being  found  on  cells  with  the  highest  PD-1  levels  (Figure  9D).  Thus,  CD25+FoxP3- 
TIL  exhibited  an  exhausted  phenotype  based  on  the  expression  of  these  markers  and 
deficient  cytokine  production. 
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To  investigate  the  prognostic  significance  of  TIL  subsets,  a  187-case  HGSC  TMA 
was  stained  with  antibodies  to  CDS,  CD25  and  FoxP3  (Figure  10).  Because  CD4  is  also 
expressed  by  macrophages,  it  is  difficult  to  score  CD4+  TIL  directly.  Instead,  we  assumed 
that  any  CD25+  and^r  FoxP3+  cell  that  did  not  express  CDS  was  a  CD4+  T  cell,  an 
assumption  that  was  justified  by  flow  cytometry  (not  shown).  With  this  staining  combination, 
we  could  directly  visualize  all  CD8+  TIL  subsets,  as  well  as  the  CD25+FoxP3-, 
CD25+FoxP3+  and  CD25-FoxP3+  subsets  of  CD4+  TIL.  To  infer  the  number  of  CD4+CD25- 
FoxP3-  TIL,  we  employed  a  second  IHC  combination  involving  antibodies  to  CDS,  CDS  and 
FoxP3.  We  determined  the  number  of  CD4+FoxP3-  cells,  which  appeared  as  CD3+CD8- 
FoxPS-  cells.  From  this,  we  subtracted  the  number  of  CD25+FoxP3-  cells  determined  from 
the  first  staining  combination.  This  yielded  an  estimate  of  the  number  of  CD4+CD25-FoxP3- 
cells.  Our  detection  and  scoring  approach  proved  valid,  as  the  results  obtained  by  multi¬ 
color  IHC  and  flow  cytometry  were  concordant.  We  stained  an  adjacent  section  of  the  TMA 
for  cytokeratin  to  unequivocally  identify  tumor  epithelium  versus  stroma;  for  all  subsequent 
analyses  we  focused  on  intraepithelial  TIL,  as  these  have  the  greatest  prognostic 
significance. 

As  expected,  CD8+  TIL  were  strongly  associated  with  disease-specific  survival  (HR 
=  0.33,  P-value<0.0001;  not  shown)  and  progression  free  survival  (HR  =  0.38,  P- 
value<0.0001)  (Figure  10A).  The  strong  prognostic  effect  of  CD8+  TIL  can  confound  the 
analysis  of  closely  associated  TIL  subsets,  such  as  Tregs.  Therefore,  to  determine  the 
prognostic  significance  of  CD4+  TIL  subsets,  we  restricted  our  analysis  to  cases  that  were 
positive  for  CD8+  TIL.  As  expected,  cases  with  a  higher  than  median  ratio  of  Tregs  to  CD8+ 
TIL  trended  toward  decreased  survival  (HR  =  1.55,  P-value=0.09)  (Figure  10B).  In  contrast, 
a  higher  than  median  ratio  of  CD25+FoxP3-  TIL  to  CD8+  TIL  was  strongly  associated  with 
survival  (HR  =  0.56,  P-value=0.02)  (Figure  IOC).  The  prognostic  effect  of  CD25+FoxP3-  TIL 
was  even  more  pronounced  in  cases  with  lower  levels  of  Tregs  (HR  =  0.29,  P-value=0.003) 
(Figure  10D).  No  other  CD4+  TIL  subsets  showed  an  association  with  survival  (not  shown). 
Thus,  the  CD25+FoxP3-  subset  was  unique  among  CD4+  TIL  in  showing  a  positive 
association  with  patient  survival. 

Thus,  we  have  discovered  a  novel  subset  of  CD4+  TIL  that,  like  FoxP3+  TIL,  express 
the  activation  marker  CD25.  Importantly,  however,  whereas  CD25+FoxP3+  TIL  are 
associated  with  poor  prognosis  (consistent  with  them  being  Tregs),  CD25+FoxP3-  T  cells 
are  associated  with  favorable  prognosis.  We  now  believe  that  these  cells  are  responsible  for 
the  favorable  prognosis  that  we  previously  attributed  to  FoxP3+  TIL.  Thus,  we  intend  to 
focus  our  antigen  discovery  efforts  on  this  subset  of  CD4+  TIL  using  the  novel  high 
througput  screen  described  in  Task  3. 

A  manuscript  describing  the  above  findings  was  recently  published  in  the  AACR 
journal  Cancer  Immunology  Research. 


KEY  RESEARCH  ACCOMPLISHMENTS: 

1)  We  published  two  manuscripts  describing  the  results  of  the  original  IDEA  award: 

(a)  Castellarin  M,  Milne  K,  Zeng  T,  Tse  K,  Mayo  M,  Zhao  Y,  Webb  JR,  Watson  PH, 
Nelson  BH,  Holt  RA.  Clonal  evolution  of  high-grade  serous  ovarian  carcinoma  from 
primary  to  recurrent  disease.  J  Pathol.  2013  Mar;229(4):515-24.  doi:1 0.1 002/path. 41 05. 
Epub  2012  Nov  29.  PubMed  PMID:  22996961. 

(b)  Wick,  D.,  Webb,  J.R.,  Nielsen,  J.S.,  Martin,  S.,  Kroeger,  D.R.,  Milne,  K.,  Castellarin, 
M.,  Twumasi-Boateng,  K.,  Watson,  P.H.,  Holt,  R.H.,  Nelson,  B.H.  2014.  Surveillance  of 
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the  tumor  mutanome  by  T  cells  during  progression  from  primary  to  recurrent  ovarian 
cancer.  Clin  Cancer  Res.  Mar  1 :20(5):1125-34.  Epub  2013  Feb  7.  PMID:  24323902. 

2)  We  published  a  key  manuscript  describing  our  discovery  of  a  prognostically  favorable  CD4+ 
T  cell  subset  in  ovarian  cancer: 

deLeeuw  RJ,  Kost  SE,  Kakal  JA,  Nelson  BH.  The  prognostic  value  of  FoxP3+  tumor- 
infiltrating  lymphocytes  in  cancer:  a  critical  review  of  the  literature.  Clin  Cancer  Res. 

2012  Jun  1:18(1 1):3022-9.  doi:  10.1 158/1078-0432.CCR-1 1-3216.  Epub  2012  Apr  17. 
Review.  PubMed  PMID:  22510350. 

3)  As  listed  below,  from  2012-2015  our  team  published  25  manuscripts  with  direct  relevance  to 
this  project,  and  58  with  indirect  relevance. 

4)  Task  1:  lllumina-based  IgG  profiling  was  completed  for  3  patients  at  3  time  points. 

5)  Task  1:  Robust  methods  were  developed  to  amplify  matched  IgG  heavy  and  light  chains  by 
single-cell  RT-PCR,  to  clone  CDR3  regions  into  IgG  expression  vectors,  and  to  prepare 
recombinant  IgG. 

6)  Task  1 :  We  discovered  that  the  bulk  of  IgG  sequences  from  tumors  are  derived  from  plasma 
cells  rather  than  CD20+  TIL.  A  manuscript  describing  these  results  is  under  review  at 
Clinical  Cancer  Research  (Appendix  A). 

7)  Task  1:  Based  on  the  above  findings,  our  IgG  cloning  and  antigen  identification  efforts  are 
now  directed  toward  these  plasma  cells  instead  of  CD20+  B  cells.  We  have  obtained  two 
new  sources  of  peer-reviewed  funding  to  continue  this  work  (see  Leveraged  Funding). 

8)  Task  2:  lllumina-based  profiling  of  TCRs  was  completed  for  3  patients  at  3  time  points. 

9)  Task  2:  Robust  methods  were  developed  to  amplify  paired  TCR-alpha  and  -beta  CDR3 
regions  by  single-cell  RT-PCR. 

10)  Task  2:  Faced  with  an  unexpectedly  high  diversity  of  TCR  sequences  among  TIL,  we 
elected  to  develop  an  entirely  new  high  throughput  method  to  identify  the  antigens 
recognized  by  TIL.  This  method  will  be  applicable  to  not  only  this  project,  but  many  future 
projects  by  us  and  others.  We  have  obtained  new  peer-reviewed  funding  to  continue  this 
work  (see  Leveraged  Funding). 

1 1 )  Task  3:  We  discovered  that  CD20+  TIL  are  highly  polyclonal,  making  it  unlikely  that  they 
serve  as  antigen  presenting  cells.  We  have  obtained  new  funding  to  investigate  the 
alternative  hypothesis  that  CD20+  TIL  serve  as  “organizers”  of  lymphocytic  infiltrates  in 
tumor  tissue,  thereby  enhancing  CD8+  TIL  responses  (see  Leveraged  Funding). 

12)  Task  4:  This  task  was  completed  in  2013/2014,  and  a  manuscript  describing  the  results  was 
published  in  Cancer  Immunology  Research  (see  Reportable  Outcomes). 


REPORTABLE  OUTCOMES: 

Manuscripts  with  direct  relevance  published  by  the  team  in  2012-2015  (n  =  25;  Teal 
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Leveraged  Funding: 

From  2012/2015,  we  obtained  six  new  grants  that  are  directiy  reievant  to  ovarian  cancer 
immunoiogy/immunotherapy: 

1.  How  does  the  immune  system  contend  with  intratumoral  heterogeneity? 

Source:  Canadian  Cancer  Society  Research  Institute 

Dates:  02/2014-01/2016  Term:  2  years 

PI:  Brad  Nelson 

Co-PIs:  Rob  Holt  and  Sohrab  Shah 

The  major  goal  of  this  project  is  to  understand  how  the  immune  system  controls  tumor 
progression  despite  the  presence  of  multiple  tumor  lineages  within  ovarian  cancer  patients. 

2.  What  proportion  of  mutations  in  the  ovarian  cancer  genome  can  be  recognized  by  the 
immune  system? 

Source:  Cancer  Research  Society 

Dates:  10/2014-09/2016  Term:  2  years 

PI:  Brad  Nelson 

The  major  goal  of  this  project  is  to  determine  whether  ovarian  tumors  express  mutated  antigens 
that  can  potentially  be  recognized  by  CD8+  T  cells  and  targeted  by  immunotherapy. 

3.  Mechanisms  of  protective  immunity  in  ovarian  cancer 
Source:  Canadian  Institutes  of  Health  Research  (CIHR) 

Dates:  10/2014-03/2015  Term:  1  year  (bridge  funding) 

PI:  Brad  Nelson 
Co-PI:  Rob  Holt 

The  major  goal  of  this  project  is  to  define  the  mechanisms  by  which  tumor-infiltrating  T  cells  and 
B  cells  work  together  to  promote  effective  tumor  immunity  in  the  setting  of  high-grade  serous 
ovarian  cancer. 

4.  Mechanisms  of  protective  immunity  in  Ovarian  Cancer 
Source:  Canadian  Institutes  of  Health  Research 

Dates:  May  2015  -  Jun  2020  Term:  5  years  (full  award) 

PI:  Brad  Nelson 
Co-PI:  Rob  Holt 

The  major  goal  of  this  project  is  to  define  the  mechanisms  by  which  tumor-infiltrating  T  cells  and 
B  cells  work  together  to  promote  effective  tumor  immunity  in  the  setting  of  high-grade  serous 
ovarian  cancer. 

5.  A  phase  I  study  evaluating  the  feasibility  and  safety  of  infusion  of  autologous  Tumour- 
Infiltrating  Lymphocytes  (TILs)  with  low-dose  interleukin-2  therapy  in  patients  with  recurrent, 
persistent  locally  advanced  metastatic  cervical  carcinoma  or  platinum-resistant  ovarian  cancer 
Source:  BC  Cancer  Foundation 

Dates:  Mar  2015  -  Feb  2018  Term:  3  years 

Co-PIs:  Brad  Nelson,  Rob  Holt  and  Anna  Tinker 
Co-applicants:  John  Webb,  Raewyn  Broday 

The  major  goal  of  this  project  is  to  conduct  a  Phase  I  clinical  trial  of  adoptive  T  cell  therapy  for 
platinum  resistant  HGSC  (and  metastatic  cervical  cancer). 
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6.  The  Immunotherapy  network  (ITNT):  Targeting  Ovarian  Cancer  (pending) 

Source:  Terry  Fox  Research  Institute 

Dates:  Mar  2016  -  Feb  2021  Term:  5  years 

PI:  Pamela  Ohashi 
Co-PI:  Brad  Nelson 

Co-applicants:  Rob  Holt  (and  12  others) 

The  major  goal  of  this  project  is  to  develop  more  effective  immunotherapies  for  HGSC  by 
defining  the  major  immune  suppressive  mechanisms  and  target  antigens  relevant  to  this 
disease. 

*Note  that  this  proposal  is  still  under  review.  However,  we  are  the  only  applicants  for  this 
funding  mechanism,  which  was  by  invitation  only.  Therefore,  there  is  virtually  a  100%  chance  of 
success. 

We  also  obtained  7  grants  with  indirect  relevance  to  this  project: 

1.  Targeting  the  Breast  Cancer  Genome  with  Personalized  Therapeutic  Vaccines 
Source:  Canadian  Breast  Cancer  Foundation  BC/Yukon 

Dates:  08/2012-07/2015  Term:  3  years 

PI:  Brad  Nelson 
Co-PI:  Rob  Holt 

The  major  goal  of  this  project  is  to  use  genome  and  transcriptome  sequencing  data  to  design 
personalized  cancer  vaccines  that  target  mutations  and  destroy  spontaneous  mammary 
tumours. 

2.  The  mutated  lymphoma  genome:  a  target  for  therapeutic  vaccination 
Source:  Canadian  Cancer  Society  Research  Institute 

Dates:  02/2013-01/2015  Term:  2  years 

PI:  Brad  Nelson 

The  major  goal  of  this  project  is  to  assess  the  immunogenicity  of  the  lymphoma  genome  and 
develop  personalized  therapeutic  vaccines. 

3.  Small-molecule  inhibitors  of  the  PD1-PDL1  interaction  for  treating  metastatic  cancer 
Source:  Canadian  Cancer  Society  Research  Institute 

Dates:  02/2013-01/2015  Term:  2  years 

PI:  Jeremy  Wulff, 

Co-applicant:  Brad  Nelson 

The  major  goal  of  this  project  is  to  develop  small  molecules  that  enhance  anti-tumor  T  cell 
responses  by  inhibiting  the  PD1-PDL1  interaction. 

4.  Mutant  MYD88:  A  target  for  adoptive  T  cell  therapy  of  WM 

Source:  The  Waldenstrom’s  Macroglobulinemia  Foundation  of  Canada  and  the  International 
Waldenstrom’s  Macroglobulinemia  Foundation  (IWMF) 

Dates:  10/2014-09/2016  Term:  2  years 

PI:  Brad  Nelson 

The  major  goal  of  this  project  is  to  develop  and  validate  TCR  contructs  that  will  allow  us  to 
engineer  CD8+  T  cells  to  specifically  recognize  a  common  mutation  in  Waldenstrom’s 
Macroglobulinemia  (MYD88L265P). 


5.  Mutation  Reactive  Autologous  T  cells  for  the  Treatment  of  Pancreatic  Cancer 
Source:  Pancreas  Centre  BC,  BC  Cancer  Foundation 
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Dates:  Sept  1, 2014  -  January  31, 2017  Term:  2.5  years 

PI:  Rob  Holt 

The  main  goal  of  this  project  is  to  develop  an  immune-based  pancreatic  ductal  adenocarcinoma 
treatment  that  will  reduce  the  chance  of  cancer  recurrence  in  these  patients. 

6.  Arming  lymphocytes  to  kill  apoptosis  resistant  cancer  cells 
Source:  Canadian  Cancer  Society  Research  Institute 

Dates:  February  1, 2015  -  January  31, 2018  Term:  3  years 

PI:  Rob  Holt 

Co-PIs:  Brad  Nelson,  Jonathan  Choy 

The  goal  of  this  project  is  to  engineer  the  granzyme-perforin  pathway  of  cytotoxic  lymphocytes 
to  deliver  to  recipient  tumor  cells  an  alternative  cytotoxic  protein  that  can  kill  by  apoptosis- 
independent  mechanisms. 

7.  Methods  and  Technology  Development  at  the  Sequencing  Platform  at  the  BC  Cancer  Agency 
Genome  Sciences  Centre 

Source:  Genome  Canada,  Genome  British  Columbia,  BC  Cancer  Foundation 
Dates:  October  1 , 2015  -  September  30,  2017  Term:  2  years 

PI:  Rob  Holt 

Co-PIs:  Marco  Marra,  Steven  Jones,  Inane  Bird,  Carl  Hansen,  Ryan  Morin  and  others 
One  of  the  goals  of  this  project  is  to  expand  the  immunogenomics  capabilities  of  the  Genome 
Canada  sequencing  platform  at  the  Genome  Sciences  Centre  by  incorporating  next-generation 
sequencing  (NGS)  based  HLA  typing  and  the  new  methods  for  NGS  based  T-cell  epitope 
discovery,  developed  as  part  of  the  current  project. 

CONCLUSION: 

Overall,  this  study  was  a  highly  successful  study.  We  discovered  two  new  subsets  of 
prognostically  important  lymphocytes  in  HGSC  (plasma  cells  and  CD25+FoxP3-  CD4+  T  cells). 
We  developed  novel,  efficient  methods  to  analyze  T  cell  and  B  cell  repertoires,  isolate  tumor- 
derived  antibodies,  and  identify  T  cell  antigens.  We  published  25  peer-reviewed  manuscripts 
with  direct  relevance  to  this  project  and  58  with  indirect  relevance.  We  competed  successfully 
for  additional  funding  to  carry  this  work  forward  and  expand  to  other  tumor  sites  (6  directly 
relevant  grants,  and  7  indirectly  relevant  grants).  Based  on  the  results  of  our  Idea  and  Teal 
Awards,  we  are  developing  a  Phase  I  clinical  trial  of  adoptive  T  cell  therapy  for  platinum- 
resistant  HGSC.  We  have  secured  over  $4M  in  funding  for  this  trial  and  the  required 
infrastructure  from  the  BC  Cancer  Foundation.  We  expect  to  treat  our  first  HGSC  patient  in  early 
2017.  We  are  very  grateful  to  the  CDMRP  for  generously  providing  both  the  Idea  and  Teal 
Awards  to  enhance  the  strength  of  our  Cancer  Immunotherapy  Program. 


APPENDICES: 

A.  Kroeger,  D.R.,  Milne,  K.,  Nelson,  BH  2015.  Plasma  cell  infiltration  identifies  ovarian 
cancers  with  robust  cytolytic  T  cell  responses  and  superior  prognosis.  Under  review  at 
Clinical  Cancer  Research. 
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None. 


FIGURES: 

See  following  pages. 
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Figure  1.  The  distribution  of  TIL  in  HGSC  and  their  association  with  long-term 

survival.  A)  Venn  diagram  showing  the  proportional  distribution  of  CD8+,  CD4+,  and 
CD20+  TIL  in  172  cases  of  optimally  de-bulked  HGSC  as  detected  by  IHC  of  a  tumor  tissue 
microarray  as  per  Milne  et  al.  2009.  B)  Kaplan-Meier  analysis  of  overall  survival  of  patients 
shown  in  panel  A.  A  small  number  of  cases  (n=5)  do  not  fit  into  the  four  categories 
described. 
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Figure  2.  Tertiary  lymphoid  structure  (TLS)  and  surrounding  TIL  in  HGSC.  Four-colour 
IHC  image  of  a  tumor  section.  CD20+  B  cells  appear  red  and  a  B  cell  follicle,  marked  “B”, 
is  found  within  the  TLS  structure.  CD3+  T  cells  appear  blue  and  the  T  cell  zone  of  the  TLS 
is  marked  “T”.  Both  B  cells  and  T  cells  are  also  found  within  the  tumor  epithelial  region, 
examples  of  which  are  marked  “TIL”.  High  endothelial  venules,  marked  “HEV”,  appear 
brown  and  were  detected  with  an  antibody  to  peripheral  node  addressin;  these  are  located 
within  the  T  cell  zone  of  the  TLS.  Activated  conventional  dendritic  cells,  “DC”,  were 
detected  with  an  antibody  to  CD208  and  appear  green-blue;  they  are  located  primarily 
within  the  T  cell  zone  of  the  TLS  (examples  are  indicated  with  arrows).  A  germinal  center, 
“GC”,  is  found  within  the  B  cell  follicle.  As  expected,  CD20  staining  intensity  is  decreased 
within  the  GC  due  to  physiological  down-regulation  of  the  CD20  protein  on  activated  GC  B 
cells. 
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Figure  3.  Tertiary  lymphoid  structures  contain  germinal  centres.  Six  color  immuno- 
histochemistry  using  antibodies  to  the  indicated  markers  was  used  to  investigate  the  cellular 
components  of  tertiary  lymphoid  structures.  We  observed  that  most  well-defined  B  cell  follicle 
structures  contained  interdigitating  networks  of  CD21+  follicular  DCs  (blue),  indicating  they 
are  germinal  centers.  Additional  staining,  not  shown  here,  demonstrated  the  presence  of 
activation-induced  cytidine  deaminase  (AID)  positive  B  cells  and  BCL-6+  T  follicular  helper 
cells  within  these  germinal  centers. 
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Figure  4.  Flow  cytometric  evidence  for  early  plasma  cells  in  HGSC  tumor 

samples.  Physically  disaggregated  human  HGSC  tumor  preparations  were  stained 
with  viability  dye,  and  then  stained  with  fluorophore-conjugated  antibodies  to  CD19, 
CD20,  CD38,  CD95  (FAS),  CD138,  IgD,  and  IgG.  After  acquisition  on  an  Influx  eight- 
channel  flow  cytometer,  data  was  analyzed  using  FloJo  software.  A)  CD19+  cells 
were  gated  from  total  viable  cells.  B)  Among  CD19+  B  cells,  we  identified  three  main 
subtypes:  naive,  memory  and  plasma  cells.  Plasma  cells  exhibited  evidence  of  class 
switching  in  that  they  were  IgD-  and  IgG''^*.  Plasma  cells  were  additionally  CD38^', 
CD20-  and  FAS+.  We  did  not  reliably  detect  CD138  expression  on  plasma  cells  by 
flow  cytometry,  similar  to  other  reports. 
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Figure  5.  Plasma  cells  arise  within  TLS  and  form  dense  stromal  infiltrates  in  HGSC. 

Serial  whole  tumor  sections  were  stained  for  the  indicated  markers.  (A)  Detection  of  TIL  and 
TLS.  TLS  appear  as  aggregates  of  CD20+  (red)  cells  containing  CD21+  follicular  DCs  (blue) 
in  germinal  centres.  (B)  Plasma  cells  appear  as  large  cells  with  CD79a+  (blue)  cytoplasm 
surrounded  by  double  positive  CD38+  (brown)  and  CD138+  (red)  membranes.  Conventional 
B  cells  express  CD79a  (blue)  on  the  cell  membrane.  CD138+  cells  in  germinal  centres 
(indicated  by  arrows)  are  thought  to  be  the  early  plasma  cells  seen  by  flow  cytometry  on  the 
preceding  figure.  Note  also  the  dense  infiltration  of  tumor  stroma  by  plasma  cells  (indicated  by 
asterisks) . 
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Figure  6.  Profiling,  cloning  and  expression  of  IqG  from  predominant  plasma  cell 

clones  in  HGSC.  (A)  Clonal  distribution  of  heavy  chain  CDR3  sequences  as  determined  by 
Sanger  sequencing  variable  portions  of  PCR-amplified  IgG  from  a  bulk  tumor  preparation 
(left),  sorted  memory  B  cells  (IgD-,  CD38-)  (centre),  or  plasma  cells  (IgD-,  CD38hi)  (right) 
from  the  same  case.  Each  pie  slice  represents  a  unique  VDJ  rearrangement,  and  the  size  of 
each  slice  reflects  the  relative  abundance  of  that  CDR3  sequence  relative  to  the  total 
population  (also  shown  by  the  numbers  below  each  pie).  Coloured  slices  indicate  CDR3 
sequences  that  are  shared  between  the  3  compartments.  Asterisks  indicate  VDJ 
rearrangements  for  which  we  have  successfully  cloned  matched  heavy  and  light  IgG  chains 
from  an  individual  plasma  cell  by  RT-PCR.  (B)  Single  tumor-infiltrating  plasma  cells  (CD19+ 
CD38hi)  were  sorted  into  individual  wells  of  a  96-well  plate.  Each  well  was  subjected  to  qRT- 
PCR  with  IgG  heavy  and  light  chain-specific  primers.  Shown  are  PCR  products 
corresponding  to  immunoglobulin  G  heavy  (IgHy)  (positive),  kappa  (IgK)  (negative),  and 
lambda  (IgL)  (positive)  from  one  representative  plasma  cell.  (C)  Human  IgG-specific 
western  blot  analysis  of  non-reduced  supernatants  from  293T  cells  co-transfected  with 
paired  heavy  and  light  chain-encoding  expression  vectors  derived  from  the  same  single 
sorted  plasma  cell  (D2)  from  (B).  Recombinant  IgG  molecules  appear  as  ~150kDa  bands. 
The  protein  ladder  appears  positive  due  to  direct  conjugation  with  fluorescent  dye.  We  have 
successfully  cloned  and  made  rIgG  from  15  single  plasma  cells. 
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Figure  7.  Molecular  cloning,  sequencing,  and  alignment  of  paired  immunoglobulin 

genes.  Following  PCR  amplification,  pairs  of  immunoglobulin  genes  were  molecularly 
cloned  into  expression  vectors.  Inserted  portions  were  then  sequenced  by  Sanger 
sequencing.  Following  sequencing,  the  variable  portions  of  cloned  immunoglobulin 
sequences  were  aligned  to  the  germline  elements  via  the  IMGTV-Quest  software. 
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Figure  8.  A  new  high  throughput  platform  for  efficient  identification  of  T  cell 

antigens.  (A)  IDS  cells  transduced  with  either  the  Ova  minigene  or  Ova-scrambled 
minigene  were  combined  at  a  1:1  ratio  to  form  a  binary  mixed  target  population.  Targets 
were  co-incubated  with  OT-I  CTL  for  4h  prior  to  FACS  analysis.  Recovered  cells  were 
lysed  and  genomic  DNA  was  used  as  template  for  qPCR  assays  using  a  custom  primer/ 
probe  set  capable  of  distinguishing  between  Ova  and  Ova-scrambled  sequences.  (B) 
Mixtures  of  IDS.minOva  and  ID8.RM  were  co-incubated  with  OT-I  CTL  at  a  1:1  target/ 
effector  ratio  for  4h  prior  to  FACS  analysis.  Recovered  cells  were  lysed  and  integrated 
minigenes  were  amplified  from  genomic  DNA  using  PCR  primers  specific  for  the 
conserved  transgene  region  flanking  the  minigene  site.  Amplicons  were  sequenced  on 
the  lllumina  MiSeq  platform  with  2x250  paired-end  chemistry. 
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Figure  9.  Phenotype  and  functional  characteristics  of  CD4+  TIL  subsets.  Multi-parameter 
flow  cytometry  was  used  to  assess  expression  of  the  indicated  activation  and  differentiation 
markers  and  cytokines  for  four  CD4+  T  cell  subsets  defined  by  CD25  and  FoxPS  expression. 

(A)  Contour  plot  showing  four  CD4+  TIL  subsets  and  color  scheme  used  in  other  panels: 
CD25+FoxP3-  (red),  CD25+FoxP3+  (orange),  CD25-FoxP3-  (blue),  CD25-FoxP3+  (green). 

(B)  Expression  of  T  cell  activation  and  differentiation  markers  by  the  four  TIL  subsets.  (C) 
Cytokine  production  after  3  hours  of  stimulation  with  PMA  and  ionomycin.  (D)  Expression  of 
exhaustion  markers  (note  that  the  right  panel  shows  only  the  CD25+FoxP3-  and  CD25-FoxP3- 
subsets).  Data  is  shown  for  one  representative  case  from  a  total  of  six.  Additional  information 
can  be  found  in  deLeeuw  et  al.  Can  Imm  Res  201 5. 
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Figure  10.  Prognostic  significance  of  TIL  subsets.  Kaplan-Meier  plots  showing  the 
association  between  TIL  subsets  and  progression  free-survival  in  a  large  HGSC  cohort. 

(A)  Density  of  intraepithelial  CD8+  TIL.  (B)  Ratio  of  CD25+FoxP3+  to  CD8+  TIL.  (C)  Ratio 
of  CD25-FoxP3+  to  CD8+  TIL.  (D)  Ratio  of  CD25+FoxP3-  to  CD8+  TIL.  In  each  panel,  log- 
rank  tests  were  used  to  determine  P-values,  and  hazard  ratios  (HR)  are  shown  for  each 
analysis.  For  panels  B-D,  analyses  were  restricted  to  cases  that  were  positive  for  CD8+ 
TIL.  Additional  information  can  be  found  in  deLeeuw  et  al.  Can  Imm  Res  2015. 
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Statement  of  Translational  Relevance: 


The  role  of  B  eells  in  antieaneer  immunity  remains  eontroversial.  Our  data  direetly 
addresses  this  eontroversy  by  demonstrating  that  plasma  eells  (PCs)  are  an  integral  eomponent  of 
robust  CD8+  tumor-infiltrating  lymphoeytes  (TIL)  responses.  We  show  that  the  well-established 
prognostie  benefit  of  CD8+  TIL  is  restrieted  to  tumors  that  additionally  harbor  PCs.  Our  findings 
indieate  that,  rather  than  working  in  opposition,  the  B  eell-  and  T  eell-lineages  mount  elosely 
integrated  responses  to  human  tumors  as  refleeted  by  their  physieal  eo-loealization,  synergistie 
funetional  profiles,  and  interdependent  prognostie  signifieanee.  We  discuss  the  implications  of 
these  findings  for  the  development  of  immunotherapies  that  engage  and  fortify  intrinsic  tumor 
surveillance  mechanisms  to  achieve  more  durable  clinical  responses. 
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Abstract: 


Purpose:  CD8+  tumor-infiltrating  lymphoyctes  (TIL)  are  key  mediators  of  anti-tumor  immunity 
and  are  strongly  associated  with  survival  in  virtually  all  solid  tumors.  However,  the  prognostic 
effect  of  CD8+  TIL  is  markedly  higher  in  the  presence  of  CD20+  B  cells,  suggesting  cooperative 
interactions  between  these  lymphocyte  subsets  lead  to  more  potent  anti-tumor  immunity. 
Experimental  Design:  We  assessed  the  co-localization  patterns,  phenotypes,  and  gene 
expression  profiles  of  tumor-associated  T-  and  B-lineage  cells  in  high  grade  serous  ovarian 
cancer  (HGSC)  by  multicolor  immunohistochemistry,  flow  cytometry,  and  bioinformatic 
analysis  of  gene  expression  data  from  The  Cancer  Genome  Atlas.  Results:  T  cells  and  B  cells 
co-localized  in  four  types  of  lymphoid  aggregate,  ranging  from  small,  diffuse  clusters  to  large, 
well-organized  tertiary  lymphoid  structures  (TLS)  resembling  activated  lymph  nodes.  TLS  were 
frequently  surrounded  by  dense  infdtrates  of  plasma  cells  (PCs),  which  comprised  up  to  90%  of 
tumor  stroma.  PCs  expressed  mature,  oligoclonal  immunoglobulin  G  transcripts,  indicative  of 
antigen-specific  responses.  PCs  were  associated  with  the  highest  levels  of  CD8+,  CD4+  and 
CD20+  TIL,  as  well  as  numerous  cytotoxicity-related  gene  products.  CD8+  TIL  carried 
prognostic  benefit  only  in  the  presence  of  PCs  and  these  other  TIL  subsets.  PCs  were 
independent  of  mutation  load,  BRCAl/2  status,  and  differentiation  antigens  but  positively 
associated  with  cancer-testis  antigens.  Conclusions:  PCs  are  associated  with  the  most  robust, 
prognostically  favorable  CD8+  TIL  responses  in  HGSC.  We  propose  that  TLS  facilitate 
coordinated  anti-tumor  responses  involving  the  combined  actions  of  cytolytic  T  cells  and 
antibody-producing  PCs. 
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Introduction: 


To  mediate  effeetive  tumor  eontrol,  the  immune  system  must  eontend  with  the  spatial 
heterogeneity  and  dynamie  evolutionary  proeesses  that  eharacterize  human  eaneer.  Advaneed 
eareinomas  are  eomposed  of  multiple  subelones  that,  despite  a  eommon  eellular  origin  and 
shared  founder  mutations,  show  divergent  genetic,  biological,  clinical,  and  immunological 
properties^  Moreover,  the  subclonal  structure  of  cancers  continually  evolves  in  response  to  the 
selective  pressures  imposed  by  the  host  and  the  cytotoxic  effects  of  treatment  .  Tumor  evolution 
can  lead  to  the  emergence  of  treatment-resistant  variants  that  ultimately  give  rise  to  fatal  disease. 
Yet  at  the  same  time,  emerging  tumor  variants  can  express  novel  antigens  that  instigate  new 
cycles  of  immune  recognition  and  attack.  Thus,  for  the  immune  system  to  successfully  contend 
with  advanced  cancers,  it  must  deploy  surveillance  and  effector  mechanisms  that  continually 
adapt  to  the  evolving  tumor.  Understanding  these  mechanisms  is  essential  for  the  development  of 
immunotherapies  that  yield  durable  responses. 

These  concepts  are  well  illustrated  by  high-grade  serous  ovarian  cancer  (HGSC)  .  At  the 
genomic  level,  HGSC  is  characterized  by  near  universal  mutations  in  the  tumor  suppressor  TP53, 
as  well  as  frequent  disruption  of  BRCAl,  BRCA2  or  other  genes  involved  in  homologous  DNA 
repair,  resulting  in  a  high  degree  of  genomic  instability"^.  Although  HGSC  has  an  intermediate 
mutation  load,  tumors  show  extensive  copy  number  variation"^.  Compounding  this  genomic 
complexity,  HGSC  disseminates  early  and  extensively  and  typically  is  not  detected  until 
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advanced  stages  .  As  a  result,  HGSC  exhibits  a  high  degree  of  spatial  heterogeneity  involving 
potentially  dozens  of  genomically  distinct  subclones  with  extensive  dissemination  across 
metastatic  sites^'^.  Although  HGSC  is  highly  sensitive  to  primary  platinum-based  chemotherapy. 
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the  development  of  ehemo-resistant  disease  is  common  and  can  bring  profound  changes  in 

6  8 

subclonal  architecture  and  mutation  profiles  '  . 

Despite  this  complexity,  there  is  a  strong  link  between  anti-tumor  immunity  and  patient 
survival  in  HGSC,  suggesting  the  immune  system  can  contend  with  tumor  heterogeneity  in  a 
substantial  proportion  of  cases.  In  particular,  the  presence  of  CD8+  tumor-infiltrating 
lymphocytes  (TIL)  in  primary  tumors  carries  a  >2-fold  increased  likelihood  of  survival^. 
Importantly,  however,  CD8+  TIL  do  not  operate  in  isolation.  We  have  shown  that  tumors 
containing  CD8+  TIL  are  often  additionally  infiltrated  by  CD20+  B  cells^*^’^^  CD20+  TIL  exhibit 
an  antigen-experienced,  immunoglobulin  G  (IgG)  positive  memory  phenotype''.  Importantly, 
tumors  containing  both  CD8+  and  CD20+  TIL  are  associated  with  higher  survival  rates  than 
those  containing  CD8+  TIL  alone.  Similar  results  have  been  reported  in  a  variety  of  other 
cancers'^"'"',  suggesting  that  effective  tumor  immunity  involves  cooperative  interactions  between 
T  cells  and  B  cells. 

In  our  prior  study",  we  found  that  CD8+  and  CD20+  TIL  often  co-localized  in  lymphoid 
aggregates  of  various  sizes  and  morphology.  This  is  reminiscent  of  autoimmune  conditions, 
where  lymphoid  aggregates  develop  in  affected  tissues.  In  rheumatoid  arthritis,  lymphoid 
aggregates  have  been  classified  into  three  grades  ranging  from  small  perivascular  collections  of 
B  and  T  cells  (grade  1)  to  large,  highly  organized  structures  resembling  lymph  nodes  (grade  3)'^. 
These  latter  aggregates,  referred  to  as  tertiary  lymphoid  structures  (TLS),  are  found  not  only  in 
autoimmunity,  but  also  in  chronic  infection,  graft  rejection,  and  cancer'^.  Like  conventional 
lymph  nodes,  TLS  harbor  prominent  B  cell  follicles  adjoined  by  discrete  T  cell  zones  containing 
CD4+  and  CD8+  T  cells,  DC,  and  high  endothelial  venules  (HEV)'^.  The  B  cell  follicles  of 
active  TLS  contain  germinal  centers  (GC)  with  interdigitating  networks  of  follicular  dendritic 
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cells  (fDC).  In  the  setting  of  eancer,  TLS  are  receiving  inereased  attention,  as  they  have  been 
assoeiated  with  favorable  prognosis  in  several  solid  tumors^^"'^. 


To  better  understand  the  mechanisms  by  which  T  cells  and  B  cells  work  together  to 
mediate  anti-tumor  immunity,  we  investigated  the  co-loealization  patterns,  phenotypes,  and  gene 
expression  profiles  of  tumor-assoeiated  T-  and  B-lineage  cells  in  HGSC.  We  found  that  the  most 
robust,  prognostically  favorable  CD8+  TIL  responses  are  aceompanied  not  only  by  CD20+  TIL 
but  by  dense  stromal  infiltrates  of  IgG+  plasma  cells.  We  propose  that  optimal  anti-tumor 
immunity  may  involve  elosely  integrated  cytolytic-  and  antibody-mediated  effector  mechanisms. 

Materials  and  Methods: 

Additional  information  is  provided  in  Supplementary  Materials. 

Patient  cohorts 

The  study  protoeol  was  approved  by  the  Research  Ethics  Board  of  the  BC  Cancer  Agency 
(BCCA)  and  University  of  British  Columbia  (Vancouver  BC).  HGSC  tumor  speeimens  were 
obtained  from  previously  untreated  patients  from  a  prospective  eohort  (treated  between  2007- 
present)  or  retrospective  cohort  (optimally  de-bulked  cases  treated  between  1984-2000;  Table 
S4)  Bioinformatic  analyses  utilized  the  HGSC  dataset  from  The  Cancer  Genome  Atlas 

Immunohistochemistry  (IHC)  and  image  analysis 

Antibodies  are  listed  in  Table  S6.  Multicolor  IHC  of  formalin-fixed  paraffin-embedded 
(FFPE)  tissue  was  performed  using  previously  described  methods^  ^  Slides  were  seanned  with  an 
Aperio  SeanScope  (Feica  Biosystems,  Wetzlar  Germany)  and  analyzed  using  ImageScope 
software  vl2.1  (Aperio  Teehnologies,  Vista  CA)  with  the  Stereology  Toolkit  v4.2.0  (ADCIS, 
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Saint-Contest,  France).  TIL  were  enumerated  in  ten  random  20x  fields,  and  cell  counts  were 
normalized  to  the  area  of  tumor  epithelium  evaluated.  PC  density  was  measured  using  a  four- 
point  scale  ;  for  survival  analyses,  cases  with  PC  scores  >1  were  deemed  positive. 

Flow  cytometry 

Disaggregated  tumor  cell  suspensions  were  washed  and  labeled  with  fluorophore- 
conjugated  monoclonal  antibodies  (Table  S6).  Flow  cytometry  and  sorting  were  performed  using 
a  BD  Influx  instrument  (Mississauga  ON). 

IgG  sequence  analysis 

RNA  was  extracted  from  FACS-purified  memory  B  cells  and  plasma  cells  (1-5  xlO 
each)  and  bulk  tumor  cells  (2.5  x  10^).  PCR  reactions  were  performed  using  primers  designed  to 
amplify  all  IGHG  variable  regions.  Up  to  192  clones  per  sample  were  subjected  to  Sanger 
sequencing. 

NanoString  gene  expression  analysis 

Total  RNA  was  prepared  from  FFPE  whole  tumor  sections  using  the  AllPrep  DNA/RNA 
FFPE  kit  (Ambion,  Eife  Technologies,  Carlsbad  CA).  Total  RNA  (200  ng)  was  analyzed  using 
the  Pan  Cancer  Immune  Profiling  panel  and  nCounter  platform  (NanoString  Technologies, 

Seattle  WA).  Data  was  normalized  using  nSolver  software. 

Bioinformatic  analysis 

The  HGSC  gene  expression  array  dataset  from  The  Cancer  Genome  Atlas  (TCGA)  was 
downloaded  from  bioconductor.org.  Corresponding  RNA-seq  data,  non-synonymous  point 
mutation  counts,  BRCAl,  BRCA2,  and  TP53  mutation  status,  and  BRCAl  promotor  methylation 
calls  were  appended. 
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Statistical  analysis 

Statistical  analyses  were  performed  using  R  v3.1.1  and  GraphPad  Prism  v6.0.  R-values  of 
less  than  0.05  were  eonsidered  signifieant  unless  otherwise  stated. 
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Results 


Multicolor  IHC  reveals  four  types  of  lymphoid  aggregate  in  HGSC 

To  visualize  the  different  lymphoid  aggregates  present  in  HGSC,  we  developed  a  6-eolor 
IHC  panel  that  enabled  simultaneous  detection  of  CD8+  and  CD4+  (CD3+CD8-)  T  cells,  CD20+ 
B  cells,  CD21+  follicular  DC  (fDC),  CD208+  activated  conventional  DC  (DC),  and  PNAd+ 
HEV-like  vessels.  The  panel  was  applied  to  whole  sections  from  30  randomly  selected  HGSC 
tumors.  To  control  for  anatomical  location,  we  stained  an  equal  number  of  specimens  from  ovary 
or  omentum;  for  7  cases,  we  stained  matched  ovary  and  omentum  samples  to  allow  direct 
comparison  between  these  sites. 

We  observed  a  variety  of  lymphoid  aggregates,  which  we  classified  into  four  types  based 
on  size,  cellular  composition,  and  degree  of  organization.  Type  1  aggregates  were  small 
(approximately  20  to  50  cells),  compact,  and  composed  of  CD4+  and  CD8+  T  cells,  B  cells,  and 
occasional  DC  (Fig,  lA),  thereby  resembling  grade  1  aggregates  in  RA^^.  Type  2  aggregates 
were  larger  (100  to  >1000  cells)  and  composed  of  CD4+  and  CD8+  T  cells  and  CD20+  B  cells 
(Fig.  IB).  These  aggregates  were  diffuse  and  lacked  discrete  zones  or  follicles.  Type  3 
aggregates  represented  fully  developed  TLS,  similar  to  grade  3  aggregates  in  RA'^.  TLS  had 
prominent  B  cell  follicles  with  germinal  center  (GC)-like  structures  characterized  by 
interdigitating  networks  of  CD21+  fDC.  In  addition,  they  contained  discrete  T  cell  zones  with 
CD4+  and  CD8+  T  cells,  DC,  and  PNAd+  HEV-like  vessels  (Fig.  1C).  Finally,  Type  4 
aggregates  were  composed  of  approximately  100-300  CD20+  B  cells  and  fDC  organized  into 
follicles  with  few  CD4-I-  and  CD8+  T  cells  (Fig.  ID).  As  these  structures  did  not  contain  clear  T 
cell  zones  and  were  primarily  (6/7)  found  in  omental  samples,  we  speculated  many  or  all  were 
normal  milky  spots  . 
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The  four  types  of  aggregate  were  strongly  assoeiated  with  one  another.  Type  1  and  2 
aggregates  were  found  in  the  same  17/30  tumors,  and  indeed  might  represent  a  eontinuum  of 
struetures.  TLS  were  found  in  7/30  tumors,  all  of  whieh  eontained  Type  1  and  2  aggregates. 

Type  4  aggregates  were  observed  in  7/30  tumors,  of  whieh  7/7  eontained  Type  1  and  2 
aggregates  and  5/7  contained  TLS.  Type  1  and  2  aggregates  were  more  prevalent  in  omental 
(1 1/23)  versus  ovarian  (6/23)  samples.  Likewise,  TLS  were  present  in  5/23  omental  versus  2/23 
ovarian  samples.  However,  none  of  these  differences  reached  statistical  significance. 

All  four  types  of  aggregate  were  strongly  associated  with  TIL.  Tumors  containing  Type  1 
or  2  aggregates  were  positive  for  CD4+  and  CD8+  TIL  in  17/17  cases  and  CD20+  TIL  in  14/17 
cases.  Likewise,  all  tumors  containing  TLS  (7/7)  contained  CD4+,  CD8+,  and  CD20+  TIL  (Fig, 
IE).  We  compared  the  density  of  TIL  in  TLS-proximal  (<500  pm  from  the  TLS  center)  versus 
TLS-distal  (>500  pm)  tumor  regions.  CD20+  TIL  showed  a  >I0-fold  greater  density  in  TLS- 
proximal  epithelium  (28.5  vs.  2.5  cells/20x  field;  P=  0.013,  Mann- Whitney  test;  Fig.  IF).  CD4-I- 
and  CD8+  TIL  showed  a  similar  trend  (72  vs.  24  and  227  vs.  59  cells/20x  field,  respectively), 
although  this  did  not  reach  statistical  significance  (Fig,  IF).  TIL  densities  were  similar  between 
omental  and  ovarian  sites  (all  P  >  0.4,  Mann- Whitney  test). 


TLS  are  associated  with  dense  plasma  cell  infiltrates  in  tumor  stroma 

To  investigate  whether  the  GC-like  structures  observed  in  TLS  were  sites  of  B  cell 
differentiation,  sections  containing  TLS  were  subjected  to  three-color  IHC  with  antibodies  to 
BCL-6,  CD3  and  CD20.  We  detected  nuclear  expression  of  BCL-6  in  both  CD20+  and  CD3+ 
cells  within  B  cell  follicles,  indicating  the  presence  of  GC  B  cells  and  Tfh,  respectively  (Fig. 

2A).  Moreover,  virtually  all  TLS  contained  AID+  CD20+  B  cells,  which  are  indicative  of 
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ongoing  immunoglobulin  class  switching  and  somatic  hypermutation  (Fig.  2B).  BCL-6+  and 
AID+  B  cells  were  also  observed  in  Type  4  aggregates  but  not  Type  1  or  2  aggregates.  Thus, 
TLS  exhibited  the  hallmarks  of  ongoing  GC  reactions. 

Tumors  were  stained  with  antibodies  to  CDS 8,  CD  138  and  CD79a  to  distinguish  PCs 
(whieh  are  CD20-  but  co-express  CD38,  CD138,  and  cytosolic  CD79a)  from  naive  and  memory 
B  eells  (whieh  are  CD38-CD138-  but  express  membranous  CD79a)  and  CD138+  tumor  or 
stromal  eells.  PCs  were  found  in  21/30  of  tumors  and  generally  formed  dense  stromal  infiltrates, 
constituting  50-90%  of  stromal  cells  (Figs.  2C,  D).  PCs  were  typieally  concentrated  near  the 
periphery  of  TLS  (Fig.  2E).  In  several  cases,  CD138+  eells  were  observed  in  GCs  (Fig.  2F). 
There  was  a  strong  assoeiation  between  TLS  and  the  density  of  stromal  PCs  (mean  PC  score  in 
the  presence  of  TLS  =  2.6  vs.  1.0  in  the  absence  of  TLS,  P  =  0.0004,  Mann-Whitney  test;  Fig. 
2G),  although  2/30  tumors  contained  abundant  PCs  (PC  seore  of  3)  in  the  absenee  of  TLS  (Fig. 
2G).  Thus,  TLS-associated  GCs  eould  serve  as  sites  of  PC  differentiation.  The  density  of  PCs 
was  also  positively  associated  with  the  intraepithelial  density  of  CD4-I-,  CD8+,  and  CD20+  TIL 
(all  Spearman’s  r  >0.6,  P  <  0.001;  Fig.  2F). 


Tumor-associated  PCs  are  IgG+CXCR3+  and  clonally  expanded 

The  phenotype  of  PCs  was  further  defined  by  flow  eytometry  of  disaggregated  viable 
tumor  samples  from  12  of  the  above  cases.  The  proportion  of  CD  19+  B  cells  varied  widely 
between  tumors  (mean  =  6.1%  of  viable  lymphocytes;  range  =  0.4-30).  Based  on  the  surface 
expression  of  IgD  and  CD38,  we  identified  three  major  CD19+  lymphoeyte  subsets:  IgD+IgG- 
CD38-  naive  B  eells,  IgD-IgG+CD38-  memory  B  cells,  and  IgD-IgG+CD38+  PCs  (Fig.  3 A).  As 
shown  previously^  the  majority  of  memory  B  eells  were  IgG  positive  (Fig.  3B).  PCs  comprised 
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from  5-88%  of  CD  19+  cells  and  also  expressed  surfaee  IgG  (Fig,  3B).  As  in  the  IHC  analysis, 
there  was  a  trend  towards  larger  proportions  of  PCs  in  tumors  eontaining  TLS  (mean  proportion 
of  CD19+  eells:  39%  vs.  25%),  although  2/12  tumors  had  abundant  PCs  (>30%  of  CD19+  eells) 
despite  lacking  TLS. 

CD  19+  cells  were  further  assessed  for  expression  of  B  eell  differentiation  markers.  As 
expeeted,  naive  B  eells  were  CD20+CD27-CD95-CD138-,  and  memory  B  eells  were 
CD20+CD27+CD95+CD138-  (Fig.  3C).  Also  as  expeeted,  PCs  were  CD20-CD27+CD95+  and 
showed  low  but  eonsistent  expression  of  CD138  (Fig,  3C).  The  chemokine  receptors  CXCR5 
and  CXCR3  reeruit  B-lineage  eells  into  follieles  and  sites  of  inflammation,  respeetively.  As 
expeeted,  naive  B  cells  expressed  CXCR5  but  not  CXCR3,  whereas  the  memory  B  eell 
population  was  a  mixture  of  CXCR5+  and  CXCR3+  eells  (Fig,  3C).  In  eontrast,  PCs  universally 
expressed  CXCR3  but  not  CXCR5.  Thus,  PCs  appeared  to  be  at  an  early  stage  of  differentiation 
and  had  an  inflammatory  ehemokine  reeeptor  profile. 

To  assess  elonality,  we  FACS-purified  PCs  from  3  tumors  and  sequeneed  the  Ig  heavy- 
ehainvariable  regions.  Bulk  tumor  samples  and  FACS-purified  memory  B  eells  were  analyzed 
for  eomparison.  In  3/3  tumors,  PC-derived  sequenees  were  restrieted  to  10-28  distinct  VDJ 
families  (Fig.  3D,  Table  SI),  suggesting  a  high  degree  of  elonal  expansion.  There  was  also 
widespread  evidenee  of  somatie  hypermutation  within  FDJ  families,  with  most  sequenees 
showing  at  least  5%  divergence  from  germline.  Although  CDR3  sequenees  from  PCs  and 
memory  B  eells  showed  some  overlap,  PC-derived  sequenees  were  most  similar  to  those  from 
bulk  tumor  in  terms  of  both  diversity  and  prevalenee  (Fig.  3D).  Thus,  PCs  were  elonally 
expanded  and  represented  the  predominant  source  of  IgG  mRNA  in  tumors. 
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A  coordinated  antitumor  response  including  PCs  is  associated  with  survival  in  HGSC 

To  assess  the  prognostic  significance  of  tumor-associated  PCs,  we  analyzed  a 
retrospective  HGSC  tissue  microarray  which  has  been  previously  evaluated  for  numerous  TIL 
subsets  ’  ’  .  PCs  were  detected  with  antibodies  to  CD138  and  CD79a,  as  our  previous  analysis 
showed  that  CD79a+CD138+  PCs  universally  expressed  CD38  (Fig,  2D,  3A).  Consistent  with 
our  initial  IHC  analysis,  PCs  were  found  in  36%  (62/172)  of  cases  and  showed  a  predominantly 
stromal  location.  When  PCs  were  considered  in  relation  to  other  TIL  subsets,  90%  of  cases 
(I55/I72)  fell  into  one  of  six  subgroups:  (a)  no  TIL  (11%  of  cases);  (b)  CD8+  TIL  alone  (7.5%); 
(c)  CD8+  and  CD4+  TIL  (20%);  (d)  CD8+  and  CD4+  TIL  with  PCs  (13%);  (e)  CD8+,  CD4+ 
and  CD20+  TIL  (16%);  and  (f)  CD8+,  CD4+  and  CD20+  TIL  with  PCs  (23%)  (Fig,  4A).  These 
subgroups  showed  no  association  with  patient  age  or  stage  of  disease  (ANOVA  P=0.777,  chi- 
squared  P=0.6I73,  respectively). 

The  six  subgroups  carried  distinct  prognostic  significance  (Fig,  4B).  The  survival  rate 
associated  with  CD8+  TIL  alone  was  similar  to  tumors  lacking  TIL  altogether.  Of  tumors 
containing  CD8+  TIL,  those  that  additionally  contained  CD4+  TIL,  CD20+  TIL  or  PCs  were 
associated  with  minor,  statistically  insignificant  increases  in  survival.  However,  tumors 
containing  CD8+,  CD4+  and  CD20+  TIL  together  with  PCs  were  associated  with  markedly 
increased  survival,  with  approximately  65%  of  patients  alive  at  10  years.  Thus,  the  prognostic 
benefit  of  TIL  in  HGSC  is  associated  with  the  combination  of  these  four  lymphocyte  subsets. 


Identification  of  a  minimal  plasma  cell  gene  expression  signature 
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To  investigate  the  relationship  between  PCs  and  the  underlying  molecular  and  genetic 
features  of  HGSC,  we  sought  to  identify  a  PC-associated  gene  expression  signature  that  could  be 
used  to  interrogate  the  TCGA  dataset.  We  assessed  the  expression  of  770  cancer/immune-related 
genes  in  19  HGSC  samples  that,  in  the  above  IHC  analyses,  were  found  to  be  (a)  positive  for 
both  PCs  and  CD20+  B  cells  (n=10),  (b)  positive  for  PCs  but  negative  for  CD20+  B  cells  (n=4), 
(c)  negative  for  PCs  but  positive  for  CD20+  B  cells  (n=2),  and  (d)  negative  for  both  PCs  and 
CD20+  B  cells  (n=3).  All  tumors  contained  CD4-I-  and  CD8+  TIL.  Average  transcript  counts 
between  the  4  subgroups  were  highly  concordant  (Fig  5A;  Pearson  r  >  0.85).  However,  66  genes 
were  expressed  at  >10-fold  higher  levels  in  PC  high  tumors  (PC  score  of  3;  n=4)  (Table  S2), 
including  the  B-lineage  genes  CD79A,  MS4A1  (CD20)  and  TNFRSF17  (B  cell  maturation 
antigen;  BCMA)  (Fig,  5A).  Of  these,  expression  of  TNFRSF17  was  uniquely  associated  with  PC 
positive  tumors  compared  to  tumors  containing  B  cells  but  not  PCs,  or  the  other  subgroups  of 
tumors  (Table  S2).  TNFRSF17/BCMA  is  known  to  be  expressed  when  B  cells  differentiate  into 
antibody-secreting  cells^"^,  and  it  is  essential  for  the  survival  of  long-lived  PCs  but  not  memory  B 
cells  .  Moreover,  TNFRSFl  7  has  been  described  as  a  signature  gene  for  human  PCs  . 

Given  our  observation  that  PCs  contribute  the  majority  of  IgG  mRNA  in  HGSC  tumors 
(Fig,  3D),  we  reasoned  that  TNFRSFl  7  and  IgG  gene  expression  levels  should  be  correlated. 
Indeed,  analysis  of  TCGA  microarray  data  revealed  a  strong  correlation  between  the  expression 
of  TNFRSFl  7  and  the  immunoglobulin  joining  region  (IGJ)  gene  segment,  which  is  common  to 
all  antibody  mRNA  transcripts  (Spearman’s  r  =  0.86;  Pearson’s  r  =  0.81).  When  we  plotted 
normalized  gene  expression  values  for  these  two  genes,  three  patient  subgroups  emerged  (Fig, 
5B):  (a)  high  expression  of  both  TNFRSFl  7  and  IGJ  (n=56),  (b)  moderate  to  high  expression  of 
IGJ hwi  not  TNFRSFl  7  (n=313),  and  (c)  negligible  expression  of  TNFRSFl  7  and  IGJ  (n=201).  A 
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similar  pattern  was  seen  using  RNA-seq  data  from  a  subset  of  these  tumors  (Fig.  5C).  As 
expected,  tumors  that  expressed  both  TNFRSFl  7  and  IGJ  also  expressed  high  levels  of  IGHV 
segments  (Fig.  5D).  The  majority  of  such  cases  expressed  5-15  IGVH  segments  (Fig.  5D), 
consistent  with  the  oligoclonal  nature  of  IgG  transcripts  we  observed  by  Sanger  sequencing  of 
IGHV  regions  (Fig.  3D).  Finally,  in  accord  with  our  flow  cytometry  data  (Fig.  3B),  IgG-derived 
transcripts  were  far  more  abundant  than  other  immunoglobulin  subtypes.  Specifically,  IGGl, 
IGG2,  and  IGG3  were  predominant,  with  low  expression  of  IGG4,  IGAl,  and  IGM  in  some  cases 
(Fig.  SI).  Thus,  TNFRSFl  7  and  IGJ  comprised  a  two-gene  signature  that  distinguished  TCGA 
cases  with  (a)  both  PCs  and  CD20+  TIL  {TNFRSFl  and/G/'^);  (b)  CD20+  TIL  without 
PCs  {TNFRSFl  and  and  (c)  neither  CD20+  TIL  nor  PCs  {TNFRSFl  and  /G/"'^). 


Association  of  the  PC  signature  with  immune-related  genes,  patient  survival,  and  CT  antigens 

The  TNFRSFl  7! IGJ  gene  signature  was  used  to  explore  associations  between  PCs  and 
other  immune-related  gene  products  in  the  TCGA  gene  expression  microarray  dataset.  Consistent 
with  our  IHC  data  (Fig.  2E),  tumors  with  a  PC  signature  (i.e.,  TNFRSFl  7^‘^^  and  IGJ’'^^)  showed 
elevated  expression  of  CXCL13  (which  contributes  to  the  establishment  and  organization  of 
TLS)  and  IL21  (which  is  produced  primarily  by  Tfh  cells)  (Fig.  6A).  Although  a  marker  of  TLS, 
CXCL13  is  unlikely  to  serve  as  a  chemoattractant  for  PCs,  as  PCs  express  the  inflammation- 
associated  chemokine  receptor  CXCR3  rather  than  CXCR5  (Fig.  3C)^^.  Indeed,  the  three 
CXCR3  ligands  {CXCL9,  CXCLIO  and  CXCLll)  were  highly  expressed  in  tumors  with  PC  or  B 
cell  signatures  (Fig.  6B),  providing  a  plausible  mechanism  for  attraction  of  PCs  to  the  tumor 
microenvironment.  Tumors  with  a  PC  signature  also  showed  high  expression  of  TNFSF13B 
(BAFF,  a  ligand  for  BCMA)  (Fig.  6C).  In  contrast,  IL6  levels  were  only  modestly  elevated,  and 
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TNFSF13  (APRIL)  and  CXCL12  levels  were  similar  to  the  PC-negative  subgroups  (Fig,  6C). 
Thus,  the  BAFF/BCMA  axis  might  provide  the  primary  growth  and  survival  signal  for  PCs  in  the 
HGSC  mieroenvironment. 

We  also  evaluated  the  relationship  between  PCs,  cytolytic  gene  signatures,  and  prognosis. 
Consistent  with  our  IHC  data  (Fig,  2H)  we  observed  a  stepwise  increase  in  the  expression  of 
CD8A  (a  marker  of  CD8+  TIL)  in  cases  with  B  cell  and  PC  signatures  (Fig,  6D),  whereas 
expression  of  the  Treg-associated  transcription  factor  FOXP3  was  uniform  across  subgroups 
(Fig,  6E).  Cases  with  the  PC  gene  signature  also  showed  elevated  expression  of  IFNG,  GZMB, 
and  PRFl  (Fig,  6F).  In  accord  with  this  cytotoxic  T  cell  signature,  the  PC  gene  signature  was 
strongly  associated  with  overall  survival  (no  B  cells  vs.  PC,  log-rank  F’=0. 0086;  B  cells  vs.  PC 
F*=0.0241),  whereas  IGJ  alone  carried  no  prognostic  benefit  (Fig,  6G).  Thus,  consistent  with  our 
IHC  data  (Fig,  4B),  PCs  were  associated  with  cytotoxic  immune  responses  and  patient  survival. 

To  gain  insight  into  potential  target  antigens  of  PCs,  we  evaluated  the  relationship 
between  the  PC  gene  signature  and  various  classes  of  tumor  antigen.  Non-synonymous  point 
mutations  can  give  rise  to  “neo-epitopes”  that  are  recognized  by  CD8+  and  CD4-I-  TIL^^. 
However,  neither  the  PC  nor  B  cell  gene  signatures  showed  an  association  with  the  total  number 
of  point  mutations  in  tumors  (chi-squared  P=0.9429;  Fig,  6H)  nor  predicted  HLA-A  associated 
neo-epitopes  (chi-squared  F’=0. 7675;  Fig  S2A).  Although  57? CH 7  impairment  has  been 
correlated  with  the  expression  of  TIL-related  genes  in  HGSC  ’  ,  the  PC  and  B  cell  gene 
signatures  showed  no  association  with  BRCAl  impairment  (including  germline  mutations, 
somatic  mutations  and  DNA  methylation)  (Chi-squared  5=0.7415;  Fig,  61)  nor  BRCA2  mutation 
(Chi-squared  5=0.9142).  Likewise,  the  PC  and  B  cell  gene  signatures  were  not  associated  with 
specific  types  of  5555  mutation  (chi-squared  5=0.7175;  Fig,  S2B).  We  also  evaluated  the 
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expression  of  45  commonly  overexpressed  antigens  and  15  well-characterized  differentiation 
antigens  (Table  S3)  but  found  no  association  with  the  PC  or  B  cell  gene  signatures  (Fig,  S2C). 
Finally,  we  assessed  the  expression  of  104  cancer-testis  (CT)  antigens.  Tumors  with  a  PC 
signature  expressed  two-fold  more  CT  antigens  on  average  than  tumors  with  a  B  cell  signature  or 
no  B  cell  signature  (cut-off:  z=3,  ANOVA  T’=0.0003;  Fig,  6J).  Similar  results  were  obtained 
using  different  thresholds  for  gene  expression  (z=2,  T’=0.0093  and  z=4,  T’<0.0001).  Semi- 
supervised  hierarchical  clustering  revealed  a  subset  of  CT  antigens  that  were  over-represented  in 
cases  with  a  PC  gene  signature  (Fig,  S3),  including  known  targets  of  autoantibody  responses 
such  as  NY-ESO-1,  MAGE  A1  and  CTAG2  (Fig,  S3;  Table  S4).  Thus,  from  this  broad 
bioinformatic  analysis,  CT  antigens  emerged  as  a  potential  class  of  target  antigen  underlying  PC 
responses  in  HGSC. 

Discussion: 

While  investigating  lymphoid  aggregates  in  HGSC,  we  discovered  a  novel, 
prognostically  significant  association  between  TIE,  TES,  and  PCs.  PCs  were  strongly  associated 
with  mature  TES  and  formed  dense  aggregates  in  tumor  stroma.  PCs  had  an  early  differentiated 
phenotype,  expressed  surface  IgG,  and  showed  evidence  of  clonal  expansion  and  somatic 
hypermutation.  PC  infiltrates  were  strongly  associated  with  CD8+  TIE  and  other  hallmarks  of 
cytotoxic  immune  responses.  Indeed,  CD8+  TIL  carried  prognostic  benefit  only  when  found  in 
combination  with  PCs,  CD20+  TIL  and  CD4-I-  TIL,  suggesting  these  four  lymphocyte  subsets 
work  in  concert  to  promote  antitumor  immunity.  Interrogation  of  the  TCGA  dataset  revealed 
correlations  between  PCs  and  other  features  of  active  antitumor  responses,  including  B-lymphoid 
growth  and  survival  factors,  TLS-associated  genes,  and  cytokines  and  chemokines  associated 
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with  cytolytic  immune  responses.  Tumors  with  a  PC  gene  signature  had  average  mutation  loads 
and  BRCAl/2  status  but  expressed  more  CT  antigens,  suggesting  the  latter  represent  target 
antigens  for  PC  responses.  Collectively,  our  findings  reveal  an  important  but  underappreciated 
collaboration  between  PCs  and  TIL  in  anti-tumor  immunity.  With  better  understanding,  it  may 
be  possible  to  enhance  these  functional  interactions  to  achieve  improved  tumor  surveillance  and 
more  durable  responses  to  immunotherapy. 

The  positive  assoeiations  we  found  between  PCs  and  tumor  immunity  are  eonsistent  with 
some  but  not  all  prior  reports  concerning  the  immunologieal  role  of  PCs.  On  the  positive  side, 
PCs  or  PC-like  gene  signatures  have  been  associated  with  favorable  prognosis  in  breast*"^,  lung^\ 
eolorectal  and  other  eaneers  .  Indeed,  in  a  reeent  pan-eancer  gene  expression  analysis,  a  PC- 
associated  gene  signature  was  among  the  strongest  positive  prognostic  factors  across  cancer 
types  .  Moreover,  prognostically  favorable  “B  cell”  gene  signatures  often  contain 
immunoglobulin  transcripts^"^’^^,  which  our  data  indicates  are  likely  attributable  to  PCs  (Fig,  3D). 

O'!  -JA 

On  the  other  hand,  both  B  eells  ’  and  PCs  have  been  shown  to  play  inhibitory  roles  in  caneer 
and  other  settings.  Recent  studies  in  infeetion  and  autoimmunity  have  identified  a  subset  of  PCs 
(“regulatory”  PCs)  that  inhibit  T  cell  responses  via  the  immunosuppressive  eytokines  IL-10  and 
IL-35  .  Shalapour  and  eolleagues  recently  deseribed  a  similar  PC  subset  in  murine  and  human 
prostate  caneer  that  inhibited  CD8+  T  eell-mediated  tumor  immunity  and  eonsequently  the 
effeetiveness  of  ehemotherapy  .  Notably,  this  PC  subset  expressed  IgA,  IL-IO  and  PD-LI  and 
differentiated  in  response  to  TGF-beta.  In  contrast,  we  found  that  PCs  in  HGSC  universally 
expressed  IgG  (Figs.  3B,  SI),  suggesting  they  undergo  class  switching  in  a  pro-immune  milieu 
dominated  by  cytokines  sueh  as  IFN-y  rather  than  TGF-p.  Moreover,  HGSC-associated  PCs 
expressed  CXCR3  (Fig.  3C),  which  is  induced  in  B-lineage  cells  by  an  IFN-y/T-bet  dependent 
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pathway  and  enables  their  reeruitment  to  inflammatory  sites  .  CXCR3+  PCs  have  also  been 

38  39 

deseribed  in  autoimmunity  ’  ,  where  they  serve  to  exaeerbate  rather  than  suppress  immune 
responses.  Thus,  the  immunosuppressive  PCs  described  in  prostate  cancer  may  be  unique  to 
that  disease,  as  the  PCs  found  in  HGSC  and  other  malignancies  have  properties  and  associations 
consistent  with  a  positive  role  in  anti-tumor  immunity. 

Given  the  strong  prognostic  significance  of  CD8+  TIL,  it  is  widely  assumed  that  they  are 
the  primary  mediators  of  anti-tumor  immunity.  However,  CD8+  TIL  are  often  functionally 
impaired  in  the  tumor  microenvironment"^'^,  raising  the  possibility  that  alternative  immune 
mechanisms  are  equally  if  not  more  important.  As  the  normal  physiological  role  of  PCs  is  to 
serve  as  “antibody  factories”,  they  could  mediate  anti-tumor  effects  by  producing  antibodies 
against  tumor-associated  antigens.  Such  effects  might  be  augmented  by  the  dense  localization  of 
PCs  in  tumor  stroma,  which  would  enable  high  concentrations  of  antibody  to  accumulate  locally. 
Theoretically,  PC-derived  antibodies  could  mediate  direct  antitumor  effects  by  binding  to  and 
disrupting  the  function  of  their  cognate  antigens,  activating  the  complement  pathway,  and/or 
triggering  antibody-dependent  cellular  cytotoxicity  (ADCC)"'\  In  this  regard,  the  predominant 
antibody  subtypes  in  HGSC  included  IgGI  and  IgG3  (Fig.  SI),  which  can  activate  both 
complement  and  ADCC"''.  Finally,  PC-derived  antibodies  could  opsonize  tumor  antigens, 
thereby  facilitating  antigen  presentation  and  broadening  of  T  ceil  responses"'^. 

Defining  the  anti-tumor  mechanisms  used  by  PCs  will  require  identification  of  their 
cognate  antigens.  Our  bioinformatic  analyses  suggest  that  PC-derived  antibodies  preferentially 
recognize  CT  antigens  over  differentiation/overexpressed  antigens  or  mutant  gene  products 
(Figs,  7,  S2).  Accordingly,  Germain  and  colleagues  recently  demonstrated  recognition  of 
multiple  CT  antigens  by  lung  cancer-associated  antibodies"'^.  Other  classes  of  antigen  might  also 
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be  relevant,  as  lung  eaneer-derived  antibodies  were  shown  to  reeognize  p53  and  various 
overexpressed  self-antigens"^"^.  Furthermore,  antibodies  derived  from  medullary  breast  eaneer 
were  shown  to  reeognize  ganglioside  D3  and  an  apoptosis-assoeiated  form  of  aetin"^^’"^^.  To  better 
understand  the  effeetor  mechanisms  used  by  PCs  in  HGSC,  we  are  cloning  IgG  molecules  from 
individual  PCs  and  attempting  to  identify  their  cognate  antigens. 

In  addition  to  antibodies,  PCs  could  influence  anti-tumor  immunity  through  cell-based 
mechanisms.  PCs  can  regulate  T  cell  responses  by  expressing  cytokines  and  other 
immunomodulatory  factors  .  Another  possibility  is  that  PCs  might  physically  exclude  immune- 
suppressive  cell  types,  such  as  cancer-associated  fibroblasts  and  myeloid-derived  suppressor 
cells"^^,  creating  a  more  permissive  tumor  microenvironment  for  CD8+  TIL  responses.  Indeed, 
lymphocyte-  and  mesenchymal-derived  gene  signatures  are  inversely  correlated  in  HCSC"^’^*^. 

Our  findings  suggest  several  novel  avenues  for  cancer  immunotherapy.  The  antibody- 
mediated  effects  of  PCs  could  potentially  be  mimicked  through  administration  of  recombinant 
antibodies  combined  with  immune  modulators  that  enhance  downstream  processes  such  as 
antigen  spreading,  CDC  and  ADCC.  The  cell-mediated  effects  of  PCs  could  potentially  be 
reproduced  by  adoptive  transfer  of  tumor-specific  PCs;  induction  of  TLS  using  vectors  that 
express  relevant  cytokines/chemokines"^^;  or  therapeutic  vaccination"^^’^*^.  By  designing 
immunotherapies  that  engage  both  the  humoral  and  cellular  arms  of  the  immune  system,  it 
should  be  possible  to  establish  endogenous  tumor  surveillance  mechanisms  that  more  effectively 
contend  with  the  spatial  heterogeneity  and  continual  evolution  of  advanced  cancers. 
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Figure  legends: 


Fig.  1.  Six-color  IHC  reveals  four  types  of  lymphoid  aggregate  in  HGSC  and  associations  with 
TIL.  Antibodies  to  CDS  (green),  CDS  (purple),  CD20  (red),  CD21  (blue),  CD208  (black),  and 
PNAd  (brown)  were  used  to  visualize  lymphoid  aggregates  in  whole  tumor  sections  from  30 
HGSC  cases.  A)  Type  1  aggregates  contained  CD8+  T  cells,  CD4+  T  cells  (detected  as 
CD3+CD8-  cells),  CD20+  B  cells,  and  CD208+  DCs.  B)  Type  2  aggregates  contained  CD8+  T 
cells,  CD4+  T  cells  and  CD20+  B  cells  in  diffuse  patterns  with  no  clear  follicles  or  T  cell  zones. 
C)  Type  3  aggregates  (TLS)  contained  B  cell  follicles  with  germinal  centers  (GC)  distinguished 
by  interdigitating  networks  of  CD21+  fDCs.  B  cell  follicles  were  adjacent  to  T  cell  zones 
containing  primarily  CD4+  T  cells  and  CD208+  DCs,  as  well  as  CD8+  T  cells.  PNAd+  vessels 
were  found  in  T  cell  zones  and  surrounding  follicles.  D)  Type  4  aggregates  contained  CD20+  B 
cells  and  CD21+  fDCs  without  clear  T  cell  zones  and  may  represent  milky  spots.  E)  Densities  of 
CD4+,  CD8+  and  CD20+  TIL  in  tumors  in  which  TLS  were  present  (n=7)  or  absent  (n=23).  F) 
Mean  densities  of  CD4+,  CD8+  and  CD20+  TIL  in  10  random  fields  versus  10  fields  within  500 
|um  of  a  TLS  center  (n=7  cases).  C-values  refer  to  means  that  were  compared  using  unpaired  (E) 
or  paired  (F)  t  tests.  Scale  bars;  A  and  D  =  50  pm,  B  and  C=100  pm. 

Fig.  2.  Evidence  of  ongoing  immune  reactions  within  TLS  in  HGSC.  TLS  were  analyzed  by 
multicolor  IHC  for  various  immune-related  markers.  A)  Three-color  stain  for  BCL-6  (brown), 
CD3  (green),  and  CD20  (red)  showing  BCL-6+  B  cells  and  follicular  helper  T  cells  (Tfh)  in  a 
GC.  B)  Two-color  stain  for  AID  (brown)  and  CD20  (red)  showing  AID+  B  cells  in  a  GC.  C) 
Three-color  stain  for  CD38  (brown),  CD79a  (blue),  and  CD138  (red)  showing  dense  stromal 
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infiltration  by  PCs  (PC  score  =  3).  D)  False-colored  image  of  PCs  showing  co-localization  of 
CD38  (green)  and  CD138  (red)  (merge  =  yellow).  E)  Example  of  a  TLS  surrounded  by  PCs.  GC 
=  germinal  center.  F)  CD138+  cells  (arrow)  within  a  TLS-associated  GC,  possibly  representing 
an  early  stage  of  PC  differentiation.  G)  Average  stromal  density  of  PCs  (on  a  seale  of  0-3)  in 
tumors  without  (n=23)  and  with  (n=7)  TLS.  E-value,  Mann- Whitney  test.  H)  Association 
between  stromal  PC  scores  and  the  average  density  of  CD4-I-,  CD8+  and  CD20+  TIL.  r  and  P 
values  refer  to  Spearmann’s  correlation.  Scale  bars:  A,  B,  D  =  50  pm;  C  =  100  pm;  E,  L  =  200 
pm. 

Fig.  3.  Cell  surface  phenotypes  and  clonality  of  B-lineage  cells  in  HGSC.  (A-C)  Multicolor  flow 
cytometry  was  used  to  analyze  TIL  from  12  disaggregated  tumor  samples.  Data  is  color-coded  as 
follows:  red  =  naive  B  cells  (IgD-l-,  CD38-),  blue  =  memory  B  cells  (IgD-  CD38-),  and  green  = 
PCs  (IgD-,  CD38hi).  A)  Representative  contour  plot  showing  expression  of  IgD  and  CD38  on 
CD  1 9+  cells.  B)  IgG  expression  on  CD  19+  cells.  C)  Expression  levels  (mean  fluorescence 
intensity,  MLI)  of  the  indicated  cell  surfaee  markers  on  CD  19+  cells  from  a  subset  tumor 
samples.  The  same  eolor-eoding  seheme  is  used  as  in  panels  A  and  B.  D)  Diversity  and  sharing 
of  Ig  V  regions  from  PCs,  memory  B  cells,  and  bulk  tumor  tissue  from  a  representative  HGSC 
ease.  The  overall  size  of  each  pie  indicates  the  relative  number  of  unique  productive  IgG 
sequences  observed  in  eaeh  sample.  The  size  of  each  pie  slice  indicates  the  relative  abundance  of 
eaeh  VDJ  rearrangement  within  a  sample.  VDJ  rearrangements  that  were  found  in  more  than  one 
sample  are  colored  to  show  the  pattern  of  sharing  between  samples.  VDJ  rearrangements  that 
were  found  in  only  one  sample  are  shown  on  a  gray  scale.  Similar  results  were  seen  with  two 
other  HGSC  cases. 
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Fig,  4.  Relationship  between  PCs,  TIL  subsets,  and  patient  survival.  A)  Venn  diagram  showing 
the  inter-relationships  between  TIL  subsets  in  the  172-case  HGSC  TMA.  Colored  circles  indicate 
patient  subgroups  that  were  positive  for  the  indicated  TIL  subsets:  blue  =  CDS,  orange  =  CD4, 
pink  =  CD20,  green  =  PCs.  Numbers  indicate  the  proportion  (%)  of  cases  in  each  subgroup.  B) 
Kaplan-Meyer  analysis  of  disease-specific  survival  for  six  patient  subgroups  based  on  the 
indicated  TIL  patterns.  C-values  refer  to  log-rank  tests  between  the  indicated  groups. 

Fig.  5.  Gene  expression  analysis  reveals  a  two-gene  PC  signature.  A)  NanoString  gene- 
expression  analysis  comparing  Logic  average  gene  counts  from  tumors  containing  no  PC 
infiltrates  (n=4)  to  those  with  heavy  PC  infiltrates  (n=3).  Highlighted  genes  are  those  expressed 
at  a  greater  than  10-fold  differential  between  groups.  B)  Standardized  microarray  gene 
expression  values  (z-scores)  for  IGJ  vs  TNFRSF17  from  TCGA  Affymetrix  U133HT  Genechip 
data  showing  three  distinct  patient  subgroups  (n=570).  C)  Standardized  read-count  data 
(fragments  per  kilobase  per  million  reads;  FPKM)  for  IGJ  and  TNFRSFI 7  from  the  TCGA 
RNA-seq  dataset  (n=273).  Data  points  are  colored  based  on  calls  from  microarray  data  (panel  B). 
D)  IGHV-XQgion  expression  (FPKM)  in  TCGA  cases  exhibiting  the  TNFRSFIIHGJ'PC 
signature,  the  IGJ'B  cell  signature,  or  no  B  cell  signature.  The  right  panel  shows  an  expanded 
view  oIIGHV  regions  from  PC  signature  positive  cases. 

Fig,6.  The  PC  gene  signature  is  associated  with  markers  of  active  humoral  and  cellular 
immunity,  patient  survival,  and  CT  antigen  expression.  As  in  Figure  5,  TCGA  cases  were 
stratified  into  three  groups  based  on  the  expression  of  IGJ  and  TNFRSFI  7:  red  =  no  B  cell 
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signature,  blue  =  B  cell  signature,  and  green  =  PC  signature.  A)  Expression  of  TLS -associated 
factors.  B)  Expression  of  chemokine  ligands  for  CXCR3.  C)  Expression  of  PC  survival  factors. 
D)  Expression  of  CD8A,  an  indicator  of  CD8+  TIL.  E)  Expression  of  FOXP3.  F)  Expression  of 
cytotoxicity-associated  gene  products.  Eor  panels  A-E,  E-values  refer  to  one-way  ANOVA  with 
Tukey post  hoc  comparison.  *P<0.05,  ***P<0.001,  ****P<0.0001.  G)  Kaplan-Meier  survival 
analysis  of  HGSC  cases  from  the  TCGA  dataset  (n  =  570).  **Log-rank  tests  between  indicated 
groups:  No  B  cells  vs.  PC  P  =  0.0086,  B  cells  vs.  PC  P=  0.0241.  H)  TCGA  RNA-seq  normalized 
gene  expression  values  for  IGJ  and  TNFRSFl  7  were  overlaid  with  calls  for  above  median  levels 
of  non-synonymous  point  mutations  (n=219).  E-value,  chi-squared  test.  I)  Analysis  was 
performed  as  in  panel  A  but  overlaid  with  BRCAl  alterations  (including  germline  and  somatic 
mutations  and  promotor  methylation)  and  BRCA2  mutations  (germline  and  somatic)  (n=316).  E- 
value,  chi-squared  test.  J)  Average  numbers  of  expressed  CT  antigens  (TCGA  microarray  gene 
expression  values  >3  standard  deviations  above  the  mean)  in  cases  with  no  B  cell  signature,  the 
B  cell  signature,  or  the  PC  signature  (n=570).  E-value,  one-way  ANOVA  with  Tukey  post  hoc 
comparisons;  **E<0.01,  ***P  <0.0001. 
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Figure  2 
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Figure  4 
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Figure  5 
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Figure  6 
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